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INTRODUCTION 


I.  SOLID  STATE  DEVICE  RESEARCH 


The  ionization  coefficients  of  electrons  and  holes  in  InP  have  been  determined  from  photo- 

multiplication  measurements  on  abrupt -junction,  low-leakage,  np+  avalanche  photodiodes. 

The  ionization  rate  of  holes  0  was  found  to  be  greater  than  that  for  electrons  a.  The 

ratio  varies  with  the  peak  electric  field  E  from  0/a  =  3.8  at  E  =  4.85  x  10^  V-cm  1  to 
r  r  ,  m  m 

0/a  =  2.7  at  E„  =  6.37  x  103  V-cm'1 
m 

The  electroabsorption  or  Franz-Keldysh  effect  has  been  measured  in  GaQ  ^In^  7^A8Q  48 

with  an  energy  gap  of  1.03  eV.  Absorption  coefficients  were  determined  by  measuring  the 
transmission  vs  reverse  bias  through  a  double -heterostructure  photodiode.  The  results  are 
in  good  agreement  with  theory  and  are  applicable  to  modulators  and  detectors  in  the  1.25-  to 
1.30-pm  band  of  interest  in  optical-fiber  transmission. 

A  subtraction  technique  has  been  developed  to  deduce  the  individual  signals  from  overlapped 
ones  in  Rutherford  backscattering  spectra.  The  method  eliminates  one  of  the  overlapping 
signals  and  yields  the  derivative  of  the  other  one,  from  which  the  two  individual  signals  are 
reconstructed.  When  the  technique  was  tested,  good  results  were  obtained  using  computer- 
simulated  backscattering  on  a  thin-film  model  which  consisted  of  five  layers  of  Au-Ag  alloys. 

An  analysis  of  a  high-speed  electrooptic  guided  wave  analog-to-digital  (A/D)  converter  has 
been  carried  out.  Design  relationships  are  given  for  the  electrooptic  chip,  the  laser  sampler, 
and  the  analog  signal  amplifier.  Based  on  the  analysis  and  available  waveguide  materials 
and  technology,  it  is  concluded  that  the  development  of  a  6 -bit  1-GS/sec  converter  should  be 
feasible. 


II.  QUANTUM  ELECTRONICS 

Frequency  coincidences  have  been  determined  between  doubled  C02 -laser  radiation  and  ab¬ 
sorption  lines  of  CO  and  NO,  and  the  anticipated  energy  requirements  have  been  calculated 
for  remote  sensing  of  these  molecules  by  topographic  reflection  using  a  differential  absorp¬ 
tion  system.  The  results  indicate  that  concentrations  typical  for  vehicle  effluent  should  be 
measurable  at  reasonable  ranges  using  the  mini-TEA  laser  as  the  primary  source. 

In  order  to  characterize  the  high-performance  limitations  of  NdP5C>14  lasers,  the  dependence 
of  the  output  power  on  duty  cycle  has  been  measured  for  a  quasi-CW  laser.  No  anomalous 
mode  or  spectral  changes  were  observed  under  different  heating  or  drive  conditions;  the  CW 
output  power  limitation  is  explained  by  a  thermally  induced  population  increase  in  the  lower 
laser  level. 

The  spectroscopy  of  terbium -yttrium  pentaphosphate  crystals  has  been  investigated  to  deter¬ 
mine  the  potential  of  these  materials  as  lasers  for  the  blue-green  spectral  region.  Radiative 
and  nonradiative  rates  have  been  measured  using  KrF*  laser  excitation. 
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A  CW  Ni:MgO  laser  has  been  operated  for  the  first  time.  Nearly  6  W  of  output  at  1.32  (im  was 
obtained  at  77  K  by  optical  pumping  with  a  1.06-|im  NdtYAG  laser.  Temperature  tuning  in 
wavelength  regions  between  1.316  and  1.409  pm  was  observed,  with  quasi -CW  operation  pos¬ 
sible  up  to  a  crystal  temperature  of  235  K. 

Experimental  studies  of  the  nonlinear  susceptibility  responsible  for  the  tripling  of  C02  radia¬ 
tion  in  cryogenic  liquids  have  been  continued.  Interference  between  the  electronic  and  the 
two-photon  resonant  vibrational  contributions  to  the  third -order  susceptibility  tensor  for 
liquid  CO-O,  mixtures  has  been  used  to  measure  the  electronic  contribution.  A  value  of 
3.9  X  10"  *  Asm/V  has  been  obtained. 

Efficient  AC  Kerr  switching  of  visible  radiation  by  infrared  radiation  has  been  demonstrated 
in  liquid  02-  These  switches  are  useful  throughout  the  infrared  spectral  region  with  ultimate 
switching  speeds  of  0.1  to  0.3  psec. 

The  applicability  of  UV  laser-induced  dissociation  to  the  processing  of  microelectronic  com¬ 
ponents  is  being  explored.  Metal  films  with  micrometer  size  features  have  been  produced 
using  a  frequency-doubled  argon  ion  laser  to  photodissociate  metal  alkyl  compounds. 

A  submillimeter  modeling  system  has  been  developed  and  applied  to  imaging  and  radar  cross- 
section  measurements  of  scaled  tactical  targets.  The  dependence  of  radar  cross  sections  on 
transmitter/receiver  polarization,  resolution  element  size,  and  target  aspect  are  being  studied 
at  frequencies  of  interest  for  a  millimeter -wave  radar  system. 


III.  MATERIALS  RESEARCH 

By  using  a  simple  analytical  model  for  GaAs  solar  cells  with  the  n+/p/p+  shallow -homojunction 
structure,  good  fits  have  been  obtained  between  computer  calculations  and  experimental  data 
for  the  external  quantum  efficiency  and  conversion  efficiency  of  cells  with  different  values  of 
n+-layer  thickness.  The  calculations  not  only  yield  values  for  material  properties  of  the  GaAs 
layers  composing  the  cells,  but  will  also  permit  the  optimization  of  cell  design  parameters. 

In  experiments  on  the  crystallization  of  amorphous  semiconductor  films  by  scanning  with  the 
slit  image  of  a  CW  laser,  several  unusual  phenomena  have  been  observed  which  include  the 
formation  of  periodic  structural  features  on  the  film  surface,  pulsations  of  the  film  tempera¬ 
ture  during  scanning,  and  runaway  crystallization  of  the  entire  film  following  first  contact 
with  the  laser  image.  By  means  of  preliminary  computer  calculations  of  the  crystallization 
front  velocity,  it  has  been  shown  that  these  phenomena  can  be  attributed  to  the  liberation  of 
heat  accompanying  the  amorphous-to-crystalline  transformation. 


IV.  MICROELECTRONICS 

Uniform  crystallographic  orientation  of  500-nm-thick  silicon  films  has  been  achieved  on 
amorphous  fused  silica  substrates  by  laser  crystallization  of  amorphous  silicon  deposited 
over  surface-relief  gratings  etched  into  the  substrates.  The  gratings  had  a  square-wave 
cross  section  with  a  3.8-^im  spatial  period  and  a  100-nm  depth.  The  <100>  directions  in  the 
silicon  were  parallel  to  the  grating  and  perpendicular  to  the  substrate  plane. 
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A  SAW/CCD  programmable  matched  filter  has  Deen  demonstrated.  This  device  is  capable  of 
correlating  an  analog  SAW  signal  of  up  to  40-MHz  bandwidth  and  a  3.5-p.sec  duration  against 
300  samples  of  an  analog  reference  previously  loaded  into  a  CCD  shift  register.  This  unique 
hybrid  structure  makes  possible  real-time  processing  of  wideband  inputs  combined  with  low- 
speed,  CCD-programmable  reference  patterns. 

A  16 -tap  CCD  transversal  filter  with  tap  weights  programmable  as  4-bit  digital  words  has 
been  fabricated  and  tested.  The  device  has  been  successfully  used  as  a  programmable 
matched  filter  for  a  chirp  signal.  The  matched-filter  output  has  a  dynamic  range  of  45  dB. 

The  spectral  quantum  efficiency  of  the  100-  x  400-element  CCD  imager  being  built  for  the 
GEODSS  (Ground  Electro-Optical  Deep  Space  Surveillance)  Program  has  been  measured.  The 
imager  is  designed  to  have  a  maximum  responsivity  to  an  air  mass  2  solar  spectrum,  which 
has  a  maximum  photon  flux  at  about  700  nm.  The  quantum  efficiency  averages  approximately 
50  percent  over  a  midband  range  from  600  to  900  nm. 

V.  ANALOG  DEVICE  TECHNOLOGY 

The  major  sources  of  spurious  signals  in  memory  correlators  have  been  identified,  and  tech¬ 
niques  for  their  suppression  have  been  developed.  The  use  of  a  split  transducer  drive  in  com¬ 
bination  with  metal  shields  under  the  ends  of  the  silicon  strips  significantly  reduces  undesir¬ 
able  bulk-wave  and  edge -detection  signals.  It  has  been  shown  that  all  the  development  goals 
for  suppression  of  spurious  signals  can  be  achieved. 

MNOS  capacitors  are  being  investigated  for  use  as  analog  nonvolatile  memory  devices.  Linear 
signal  storage  over  a  dynamic  range  of  at  least  47  dB  has  been  demonstrated.  The  loss  of 
stored  charge  proceeds  logarithmically  in  time  at  a  rate  proportional  to  the  initial  stored 
charge;  a  plot  of  flat-band  voltage  vs  the  log  of  the  storage  time  is  a  straight  line  with  a  slope 
of  approximately  0.1  V/decade  per  volt  of  initial  flat-band  voltage.  Further  investigations  are 
centered  on  optimizing  the  conduction  properties  of  the  oxide  and  nitride. 
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I.  SOLID  STATE  DEVICE  RESEARCH 


A.  IONIZATION  COEFFICIENTS  OF  ELECTRONS  AND  HOLES  IN  InP 

Recent  interest  in  avalanche  photodiodes1'2  and  IMPATT  devices5  fabricated  from  InP  has 
prompted  an  investigation  of  the  ionization  coefficients  of  electrons  (a)  and  holes  (0)  in  this 
material.  The  ionization  rates  are  known  to  have  a  strong  influence  on  both  the  noise^  and 
gain-bandwidth  product  limitations5  of  avalanche  photodiodes  as  well  as  on  the  noise  and  effi¬ 
ciency  of  IMPATTs/’ 

Unambiguous  determination  of  the  ionization  rates  from  photomultiplication  data  on  ava¬ 
lanche  photodiodes  requires  two  important  experimental  conditions.7  One  is  that  pure  electron 
and  pure  hole  injection  into  the  high  field  region  must  be  obtained  in  the  same  device.  The 
second  requirement  is  that  the  magnitude  and  variation  of  the  junction  electric  field  be  known 
accurately.  These  two  requirements  are  best  met  by  use  of  an  abrupt  p-n  junction  photodiode, 
where  incident  radiation  is  completely  absorbed  in  undepleted  material  on  either  side  of  the 

junction,  allowing  either  electrons  or  holes  to  diffuse  into  the  field  region.  This  method  has 
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been  employed  successfully  to  determine  ionization  rates  in  other  III-V  materials.  ’  Earlier 
measurements  of  the  ionization  rates  in  InP  have  been  reported1^  for  diffused  junction  devices 
where  the  photomultiplication  data  were  interpreted  under  the  invalid  assumption  that  a  =  p. 
Measurements  also  have  been  reported  11  on  InP  Schottky  barrier  diodes,  a  configuration  for 
which  the  interpretation  of  the  data  is  subject  to  several  uncertainties,12  with  the  result  p/a  «  5. 

The  avalanche  photodiodes  used  here  were  of  the  inverted  mesa  geometry  fabricated  by  etch¬ 
ing  138-p.m-dia.  mesas  in  a  wafer  consisting  of  a  5-pm-thick  n-InP  layer  (n  =  1.7  x  lO1^1  cm  5) 

,  18-3 

grown  by  LPE  on  a  (lOO)-oriented  p  -InP  (Zn)  substrate  (p  ~  2  x  10  cm  ).  Mesas  were 
etched  using  a  1:1  mixture  of  3%  Br-methanol:H3P04  at  45°C.  Ohmic  n-  and  p-type  contacts 
were  made  by  microalloying  plated  AuSn  and  evaporated  AuMg,  respectively.  The  center  portion 
of  the  substrate  was  thinned  to  within  4  pm  of  the  junction  using  a  jet  of  a  1 : 1 :6 : 1  mixture  of 
HAc:HC104:HN0j :HC1  to  facilitate  photoinjection  of  electrons.15  The  final  device  structure  is 
shown  in  the  insert  in  Fig.  1-1.  The  etched  moat  around  the  base  of  the  mesa  is  an  artifact  of 
the  etching  procedure  and  has  the  positive  attribute  of  insuring  electrical  isolation  of  the  mesa 
while  retaining  some  additional  physical  support  in  the  rest  of  the  wafer. 

Dark  currents  of  the  diodes  were  of  subnanoampere  magnitudes  up  to  voltages  very  close  to 
breakdown.  Figure  1-1  illustrates  the  I-V  characteristic  of  a  typical  diode.  For  some  of  the 
better  devices,  uniform  photoresponse  gains  (within  ±10  percent)  of  about  12  were  measured 
with  leakage  currents  as  low  as  0.35  nA.  Typically,  however,  leakage  currents  at  this  value  of 
gain  tended  to  be  about  an  order  of  magnitude  larger.  Reasonably  uniform  photoresponse  gains 
of  150  with  a  concurrent  leakage  of  120  nA  have  been  measured.  Capacitance-voltage  data  in¬ 
dicated  that  these  devices  were  abrupt  junction  diodes.  Scanned  photoresponse  measurements 
from  both  sides  of  the  wafer  using  the  6 328 -A  line  of  a  HeNe  laser  were  performed  to  ensure 
the  absence  of  microplasmas  and  edge  breakdown  in  the  devices  chosen  for  study.  The  diodes 
selected  for  photomultiplication  measurements  had  uniform  multiplication  (±10  percent)  over 
the  area  of  the  diode  for  the  range  of  voltages  used  in  the  experiment. 

Photomultiplication  measurements  were  performed  using  a  focused  beam  of  strongly  ab¬ 
sorbed  radiation  (6328  A),  incident  alternately  on  the  n  and  p+  sides  of  the  diodes,  resulting  in 
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pure  hole  and  electron  injection,  respectively.  The  photoresponse  as  a  function  of  reverse  bias 
for  each  case  was  monitored  with  a  lock-in  amplifier.  The  injected  photocurrent  was  of  the 
order  of  5  nA.  Over  a  rather  wide  range  of  increasing  reverse  bias,  below  which  multiplication 
was  evident,  the  photocurrent  rose  slowly  and  very  nearly  linearly  due  to  the  increased  collec¬ 
tion  efficiency  of  the  widening  depletion  region.  A  linear  extrapolation  of  this  low  field  photo¬ 
response  was  subtracted  out  at  the  higher  fields  to  yield  the  values  for  the  hole -initiated  multi¬ 
plication  (Mp)  and  the  electron-initiated  multiplication  (Mn>  shown  in  Fig.  1-2.  The  ionization 
rates  were  determined  as  a  function  of  peak  electric  field  (Em)  from  the  multiplication  data  by 

means  of  the  following  equations  which  apply  for  the  case  of  a  one-sided,  abrupt -junction  device 

14 

with  pure  electron  and  pure  hole  injection,  respectively. 


a(Em)=^(Em)  +  Em^  jln[^-]| 

E 

n  r 


(I-D 


(1-2) 


Em  was  calculated  using  the  applied  voltage  (V)  and  the  carrier  concentration  obtained  from 
capacitance-voltage  measurements. 

Application  of  Eqs.  (1-1)  and  (1-2)  to  the  multiplication  data  of  Fig.  1-2  yields  the  hole  and 
electron  ionization  rates,  plotted  in  Fig.  1-3  as  a  function  of  l/Em-  These  results  show  clearly 
that  the  ionization  rate  for  holes  is  greater  than  that  for  electrons  over  the  entire  range  of 


Fig.  1-3.  Ionization  rates  of  electrons 
and  holes  in  InP  as  a  function  of  re¬ 
ciprocal  peak  electric  field. 
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electric  field  strengths  considered.  For  a  peak  electric  field  strength  of  4.85  x  10®  V-cm  the 
ratio  of  hole  ionization  rate  to  electron  ionization  rate  (0/a)  is  3.8.  This  ratio  decreases  with 
increasing  electric  field  until  0/a  =  2.7  at  Em  =  6.37  x  10®  V-cm  *.  A  least-squares  fit  of 
these  data  indicates  that  the  ionization  rates  can  be  approximated  over  the  range  of  electric 
fields  used  in  this  experiment  as 

a(E)  =  1.07  x  107  exp [-4.31  X  106/E)  (1-3) 

0(E)  =  9.63  X  106  exp [-3.61  X  106/E]  .  (1-4) 

The  measurements  reported  here  provide  some  initial  values  for  the  ionization  rates  in 
(lOO)-InP,  results  which  should  be  reliable  because  of  the  unambiguous  nature  of  the  experi¬ 
mental  technique  employed.  Also,  it  has  been  demonstrated  that  uniform,  high-gain,  low- 
leakage  avalanche  photodiodes,  suitable  for  use  as  detectors  as  well  as  for  multiplication  mea¬ 
surements,  can  be  fabricated.  C.  A.  Armiento 

S.  H.  Groves 
C.  E.  Hurwitz 


B.  ELECTROABSORPTION  IN  GalnAsP 

Measurements  of  electroabsorption  (the  Franz -Keldysh  effect)  have  been  made  on  the 
quaternary  semiconductor  alloy  GaQ  24Ino  76^s0  52^0.48  has  an  enerBy  EaP  Eq  =  e'^> 

The  data  were  obtained  by  measuring  the  transmission  perpendicular  to  the  plane  of  a  reverse- 
biased  double-heterojunction  GalnAsP/lnP  photodiode.  Since  the  thin  quaternary  layer  enabled 
measurements  of  absorption  coefficient  up  to  values  of  10  cm  ,  this  represents  a  significant 
extension  of  the  earlier  work  of  Stillman  et  al.3®  on  GaAs.  The  results  are  consistent  with 
electroabsorption  theory  and  are  particularly  applicable  to  detectors  and  modulators  operating 
in  the  1.25-  to  1.30-(im  wavelength  region  most  suitable  for  optical-fiber  transmission. 

The  diode  was  a  mesa  structure,  previously  described  by  Hurwitz  and  Hsieh,*^  with  a 
p+-lnP  substrate,  a  2.5-pm  p+-n  quaternary  layer,  and  an  n+-InP  top  layer  with  small  Au-Sn 
ohmic  contacts.  The  transmission  through  the  diode  was  measured  as  a  function  of  wavelength 
and  reverse  bias,  using  a  prism  monochromator  apertured  to  produce  a  focused  spot  equal  to 
the  diode  diameter  of  0.25  mm.  The  transmission  of  the  p+-InP  substrate  was  measured  by 
moving  the  spot  away  from  the  mesa;  the  ratio  of  these  two  measurements  yielded  the  trans¬ 
mission  of  the  quaternary  material  since  the  extra  InP  in  the  mesa  structure  is  effectively 
transparent  at  the  wavelengths  of  interest  and  the  off-mesa  measurement  takes  acount  of  the 
Fresnel  reflections  and  residual  absorption  in  the  thicker  substrate.  Data  were  obtained  up 
to  reverse  voltages  of  40  V,  corresponding  to  a  maximum  electric  field  of  3.6  X  10®  V/cm  for 
the  constant  donor  density  NQ  =  1.2  x  1016  cm"3  determined  from  the  capacitance-voltage 
characteristic. 

The  transmission  is  affected  only  by  the  field  distribution  in  the  n-type  material  since  the 

space-charge  layer  in  the  p+  material  is  of  negligible  thickness  because  of  the  high  acceptor 
4  8  -3 

density  10  cm  .  The  transmission  may  be  written 


In  T  = 


a(E) dx  - 


£ 


a(0)  dx 


(1-5) 


where  the  space-charge  region  extends  from  0  to  w,  and  t.  is  the  thickness  of  the  n-type 
GalnAsP.  For  constant  donor  density,  and  consequently  a  linear  decrease  of  field  with  distance. 
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a  small  change  in  reverse  voltage  results  in  a  change  in  space-charge  width  dw  and  the  change 
in  transmission  becomes 


d(lnT)  =  -a(E)  dw  +  a(0)  dw  (1-6) 

where  E  is  the  peak  field  located  at  x  =  0.  Physically,  the  whole  field  distribution  shifts  to 
higher  x;  a  thin  layer  dw  at  maximum  field  is  added  and  an  equal  thickness  layer  at  zero  field 
is  removed.  The  final  expression  in  terms  of  the  reverse  bias  becomes 

o(E)  -  a(0)  =  (dT/dw)/T  =  (ZqNjjV/ic^)1/2  (dT/dV)/T  .  (1-7) 

Using  this  expression,  the  measured  transmission  vs  reverse  bias,  a  dielectric  constant17  of 
13.2,  and  the  energy  gap,  1.03  eV,  the  absorption  coefficient  vs  field  has  been  plotted  in  Fig.  1-4. 
Also  shown  by  the  dashed  line  is  the  zero-bias  absorption  coefficient  which  was  used  to  deter¬ 
mine  o(0).  This  zero-field  value,  determined  from  the  net  transmission  through  the  whole 
quaternary  layer,  includes  the  effects  of  absorption  in  the  p+-layer.  It  therefore  introduces 
some  uncertainty  in  the  absolute  value  of  o(E)  near  the  band  edge,  although  the  difference 
[o(E)  -  o(0)]  is  still  reliable. 

WAVELENGTH  (pm) 


Fig.  1-4.  Experimental  and  theoretical  absorption  coefficient 
vs  wavelength  and  electric  field  for  GalnAsP  and  GaAs.  The 
theory  and  GaAs  data  are  from  Ref.  15. 

The  data  in  Fig.  1-4  include  the  results  for  GaAs  (Ref.  15)  as  well  as  the  theoretical  predic¬ 
tions  from  the  same  work.  The  only  materials  parameters  which  enter  into  the  theory  of  the 
1 8 

electroabsorption  effect  are  the  effective  masses  and  the  index  of  refraction.  Since  the 
19 

masses  and  index  are  the  same  within  better  than  10  percent  for  GaAs  and  GalnAsP,  the  theory 
should  apply  to  the  latter  material  also.  The  data  in  Fig.  1-4  show  good  qualitative  agreement 
with  theory,  including  the  predicted  decrease  at  energies  above  the  bandgap.  There  is  also  con¬ 
sistency  with  the  earlier  GaAs  data  in  that  the  experimental  values  are  significantly  lower  than 
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-  0.043+V 

-0.020 

+0.003 


Fig.  1-5.  Absorption  coefficient  vs 
electric  field  for  several  values  of 
(EG-ftoj)  for  Gain  As  P  and  GaAs. 
The  data  points  at  4  and  6  x  104  V/fcm 
are  for  GaAs,  Ref.  15. 


the  theory,  or  alternatively,  a  higher  electric  field  is  required  to  produce  the  predicted  absorp¬ 
tion  curve.  Figure  1-5  is  a  plot  of  absorption  coefficient  vs  electric  field  for  various  energy 
differences  below  and  above  the  gap.  Included  are  additional  data  points  on  GalnAsP  as  well  as 
points  from  the  GaAs  results.  These  data  should  be  applicable  for  GalnAsP  alloys  over  a  wide 
range  of  energy  gaps,  provided  there  is  no  significant  change  in  effective  masses  or  index  of 
refraction.  R  H  ^ngston 


C.  SUBTRACTION  OF  SIGNAL  OVERLAPS  IN  RUTHERFORD 

BACKSCATTERING  SPECTROMETRY 

Signal -overlapping  has  been  a  serious  limitation  in  Rutherford  backscattering  (RBS)  analysis 
of  thin-film  interactions.  This  problem  was  encountered  in  our  recent  RBS  study  of  Au/lnP  con¬ 
tacts20  in  which  the  Au  and  In  signals  of  the  alloyed  samples  overlapped.  To  overcome  this 
difficulty,  we  have  developed  a  technique  in  which  the  individual  signals  can  be  obtained  from 
the  overlapped  ones  by  a  proper  subtraction  between  two  RBS  spectra  of  slightly  different  inci¬ 
dent  ion  -nergies. 

Figure  1-6  illustrates  the  main  idea  of  this  technique.  In  Fig.  1-6 (a),  f^(E)  and  gj(E)  repre¬ 
sent  the  ua^kscattered  signals  of  two  different  elements  F  and  G,  respectively,  in  the  target 
which  has  been  analyzed  with  an  incident  ion  energy  Ej.  Since  f ^ ( E )  and  gj(E)  overlap,  it  is  the 
sum  of  the  two  signals,  h.  (E)  =  f.(E)  +  g,(E),  that  is  obtained  in  the  RBS  experiment.  When  a 
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Fig.  1-6.  Illustration  of  the  main  idea  in  the  present  technique. 
Parts  (a)  and  (b)  represent  RBS  spectra  of  a  sample  taken  at 
incident  ion  energies  Ej  and  E2,  respectively,  with  E2  >  Ej. 
Part  (c)  shows  the  matching  of  the  two  spectra  in  (a)  and  (b). 


slightly  higher  incident  ion  energy  E2  is  used  for  the  same  experiment,  the  signals  will  appear 
in  higher-energy  regions  as  shown  in  Fig.  I -6(b).  We  first  assume  that  the  energy  dependence 
of  ion  energy  loss  in  the  target  material  is  so  small  that  the  individual  signals  are  shifted  with¬ 
out  changes  in  shape,  l.e., 

f2(E)  =  f^E  -KpAE)  (I-8a) 

g2(E)  =  gl(B  -KqAE)  (I— 8b) 

where  Kp,  and  K^,  are  the  backscattering  kinematic  factors*1  of  elements  F  and  G.  respectively, 
and  AE  s  E2  -  Ej.  By  defining  E  5  E  +  KpAE.  Eqs.  (I-8a)  and  (I-8b)  are  then  used  to  transform 
h2(E)  =  f2(E)  ♦  g2(E)  into 

h2(E)  »  fjfE)  ♦  gt(E  *  «E) 
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(1-9) 


where  6E  5  (K^,  -  K^)  AE.  Because  6E  is  relatively  small,  Eq.  (1-9)  can  be  further  transformed 
into 


g,(E) 


I 


h2(E)  -ht(E) 


6E 


dE 


(1-10) 


The  physical  picture  is  illustrated  in  Fig.  I-6(c).  The  function  h^fE)  actually  represents  a  trans¬ 
lation  of  h2(E)  to  a  lower-energy  region  by  a  distance  KpAE,  which  makes  f2(E)  coincide  com¬ 
pletely  with  fj(E).  But  g2(E)  is  mismatched  with  ( E )  by  a  small  energy  distance  6E.  There¬ 
fore,  the  subtraction  between  h2(S)  and  hj(E)  yields  the  derivative  of  g1(E)  from  which  gj(E)  can 
be  reconstructed. 

Equations  (I-8a)  and  (I-8b)  neglect  the  energy  dependence  of  ion  energy  loss  in  the  target 
material.  Over  a  fairly  large  energy  range  of  interest  the  spectral  width  of  an  elemental  thin- 
film  target  decreases  (with  increasing  incident  ion  energy)  at  a  constant  rate.2*’22  We  assume 
that  this  is  also  true  for  signals  of  a  compound  target.  Thus,  in  addition  to  a  shift  by  a  distance 
KpAE,  the  energy  scale  of  h2(E)  should  then  be  expanded  by  a  factor  of  (1  +  sAE)  such  that  f2(E) 
matches  with  fj(E)  completely.  (The  coefficient  s  is  the  constant  "shrink  rate"  of  the  signals.) 
Because  of  this  expansion,  the  relative  energy  shift  6E  between  gj(E)  and  g.,(E)  will  increase  by 
an  amount  (K^,  -  K^)  E2  •  sAE.  Thus,  the  total  6E  is  given  by 


<5E  =  (Kp  -Kq)  AE  +  (Kp  —  Kg)  E2  •  SAE  .  (1-11) 

23 

To  test  this  technique,  we  present  a  numerical  example  in  which  a  computer  program  is 

used  to  simulate  the  backscattering  experiment.  Backscattering  spectra  of  2.0  and  2.1  MeV 
4  +  o 

He  ions  are  obtained  for  a  sample  consisting  of  5  layers  (400  A  each)  of  AuxAgt  ,  as  shown 

in  Fig.  I-7(a).  The  4He+  dose  used  for  the  2.1-MeV  spectrum  is  1.1  times  that  of  the  2.0-MeV 


Fig.  1-7.  Example  of  the  present  technique. 
Part  (a)  shows  computer-simulated  RBS  spec¬ 
tra  of  the  AuxAgj.x  sample.  Part  (b)  shows 
the  difference  between  the  two  matched  spec¬ 
tra  in  (a).  Part  (c)  shows  the  reconstructed 
as  well  as  the  original  Ag  and  Au  signals. 


8 


i 


■ 


one  (to  compensate  the  energy  dependence  of  the  Rutherford  cross  section),  such  that  the  "clean" 

Au  signals  in  the  high-energy  regions  are  of  equal  heights.  The  energy  scale  for  the  2.1  MeV 

spectrum  is  3.92  keV/channel  as  compared  to  4.00  keV /channel  for  the  2.0-MeV  one.  This 

-4  -2 

corresponds  to  an  s-value  of  2.0  x  10  /keV  or  sAE  =  2.0  x  10  .  The  SE  as  obtained  by  using 

Eq.  (X - 1 1  >  is  8.28  keV.  Figure  I-7(b)  shows  the  difference  between  the  two  matched  spectra  in 
Fig.  I-7(a).  Starting  from  the  high-energy  region,  the  integration  of  Eq.  (I - 1 1 )  is  carried  out  by 
summing  the  readings  in  Fig.  I-7(b),  channel  by  channel.  This  integration  yields  the  recon¬ 
structed  Ag  signal  shown  in  Fig.  I  -  7  ( c ).  The  reconstructed  Au  signal  is  then  obtained  as  the  dif¬ 
ference  between  the  2.0-MeV  spectrum  in  Fig.  I-7(a)  and  the  reconstructed  Ag  signal  in  Fig.  I-7(c). 

The  computer  generates  the  individual  Ag  and  Au  signals  as  well  as  the  overlapped  ones  in 
Fig.  I-7(a).  Comparisons  in  Fig.  I  -  7  (c )  show  that  the  present  technique  has  closely  reproduced 
the  computer  generated  signals  in  large  portions  of  the  spectrum.  Larger  discrepancies  are 
found  in  the  reconstructed  Au  signal  near  the  Ag  edge.  This  is  probably  due  to  the  rapid  changes 
of  the  Ag  signal  near  the  Ag  edge  which  cause  larger  errors  in  approximating  g^(E  +  6E)  —  g1(E) 
as  (dg^/dE)  6E  in  Eq.  (I- 1 0 ).  There  is  also  a  gradual  increase  in  discrepancy  between  the  two 
Ag  signals  in  the  lower  energy  region  in  Fig.  I-7(c).  This  can  be  caused  by  the  slightly  faster 
rising  rate  of  the  Ag  signal  of  the  2.0-MeV  spectrum  than  that  of  the  2.1-MeV  one,  due  to  the 
energy  dependence  of  the  Rutherford  cross  section.  Although  the  difference  is  very  small 
[causing  ~0. 5-percent  difference  in  the  two  corresponding  Ag  signal  heights  near  the  low  energy 
edges  of  the  spectra  in  Fig.  I-7(a)],  the  error  can  accumulate  because  of  the  integration  proce¬ 
dure  used  in  obtaining  the  reconstructed  Ag  signal.  Calculations  show  that  the  error  rises 
nearly  parabolically,  with  a  magnitude  of  a  few  percent  near  the  low  energy  edge  of  the  Ag  sig¬ 
nal,  which  is  comparable  to  that  observed  in  Fig.  I-7(c). 

One  important  approximation  used  in  this  technique  is  that  the  "shrink"  of  signals  can  be 
corrected  by  a  single  parameter  s.  Although  this  approximation  appears  to  be  a  good  one  in  the 
case  of  Fig.  1-7,  we  do  not  know  how  good  it  is  in  other  systems,  especially  those  with  extreme 
compositional  nonuniformity. 

In  conclusion,  we  have  demonstrated  that  individual  signals  can  be  obtained  to  high  accura¬ 
cies  from  the  overlapped  ones,  by  using  a  simple  analytical  technique  which  is  mainly  based  on 
backscattering  kinematics.  The  general  validity  of  the  approximation  that  the  dependence  of 
signal  shapes  on  the  incident  ion  energy  can  be  compensated  by  a  simple  change  of  energy  scale 
is  yet  to  be  evaluated.  Should  this  question  be  answered  affirmatively,  the  present  technique 
can  become  very  useful  in  RBS  study  of  materials  with  moderate  mass  differences  between  the 


constituent  elements. 


Z.  L.  Liau 


D.  DESIGN  CONSIDERATIONS  FOR  A  HIGH-SPEED  ELECTROOPTIC 

ANALOG-TO-DIGITAL  CONVERTER 

One  of  the  most  promising  systems  applications  of  integrated  optical  devices  is  the  utiliza¬ 
tion  of  an  array  of  waveguide  electrooptic  modulators  to  perform  analog-to-digital  (A/D)  conver- 
sion.  Electrooptic  A/D  converters  are  of  special  interest  because  they  should  be  capable 

of  performing  conversions  at  rates  far  exceeding  what  can  presently  be  achieved  with  state-of- 
the-art  technology.  In  addition,  the  electrooptic  A/D  converter  permits  optical  sampling  of  the 
analog  waveform  and  incorporates  techniques  that  eliminate  some  of  the  inherent  problems  and 
complexities  of  integrated-circuit  A/D  converters.  In  this  section,  detailed  design  considera¬ 
tions  are  presented  which  suggest  the  feasibility  of  developing  a  6 -bit  1-GS/sec  electrooptic 
A/D  converter. 
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1.  Converter  Operational  Principles 

The  fundamental  integrated  optical  component  of  the  electrooptic  A/D  converter  is  a  planar 

27-29 

waveguide  version  of  a  Mach-Zehnder  interferometric  modulator,  as  shown  schematically 
in  Fig.  1-8.  It  consists  of  an  electrooptic  crystal  (e.  g.t  LiTaQj)  containing  a  single-mode  input 


Fig.  1-8.  Schematic  drawing  of  integrated 
optical  Mach-Zehnder  interferometer  mod¬ 
ulator  fabricated  from  single-mode  wave¬ 
guides  in  LiTaOj. 


optical  waveguide  which  branches  at  a  Y  to  s;>lit  the  optical  power  into  two  equal  components. 
Each  component  of  light  travels  an  equal  distance  before  recombining  at  a  second  Y  and  exiting 
the  crystal.  In  the  absence  of  an  applied  voltage  the  light  components  recombine  in  phase,  and 
a  maximum  output  is  obtained.  If  voltage  is  applied  to  one  arm  of  the  modulator  via  the  coplanar 
electrodes  shown  in  the  figure  the  phase  velocity  of  light  propagating  in  that  arm  will  be  altered 
due  to  the  linear  electrooptic  effect.  This  will  produce  a  phase  difference  between  light  in  the 
two  modulator  arms  when  they  recombine,  which  will  reduce  the  output  signal.  In  the  extreme 
case  the  two  waves  will  be  ir  radians  out  of  phase.  This  will  result  in  a  null  output  since  the 
two  recombining  light  waves  will  try  to  form  the  second-order  mode  which  the  single-mode  guide 
cannot  support,  and  the  light  will  radiate  into  the  substrate.  In  general,  it  can  be  shown  that  for 
an  input  intensity  I  ,  the  output  intensity  varies  as 

I  =  IQ  [1  +  cos<«p  +  *s)]/2  (1-12) 

where  <p  is  the  net  phase  difference  between  the  two  guides  due  to  the  applied  voltage  and 
accounts  for  any  static  phase  difference  between  the  arms.  It  should  also  be  noted  that  the  use 
of  a  three -electrode  configuration,  one  between  the  two  waveguides  and  two  outside  the  guides, 
will  reduce  the  required  drive  voltage  by  a  factor  of  2  for  a  given  ip  and  minimize  the  modulator 
temperature  dependence. 

One  can  take  advantage  of  the  sinusoidal  variation  of  the  output  intensity  with  voltage  to  form 
an  A/D  converter  because  the  value  of  a  bit  in  a  binary  representation  of  an  analog  signal  varies 
periodically  with  thf  value  of  that  signal.  In  particular,  by  combining  a  number  of  these  modu¬ 
lators  in  parallel  as  shown  in  Fig.  I-9(a),  a  converter  can  be  fabricated  which  only  requires  one 
modulator  per  bit  of  ,-esolution.  The  intensity  output  of  the  modulators  vs  the  applied  analog 
signal  V&  is  shown  in  Fig.  I-9(b).  The  binary  code  is  obtained  by  comparing  the  modulator  out¬ 
puts  as  detected  by  a  photodiode  at  each  output  with  a  reference  signal  from  the  same  light 
source  used  to  excite  the  modulators.  As  shown  in  the  figure,  if  the  output  intensity  is  above 
or  below  threshold  It>  a  "0"  or  "1".  respectively,  is  generated.  The  corresponding  binary 
Gray  code  is  also  shown  in  the  figurts.  In  this  converter,  the  sampling  of  the  analog  signal 
is  done  optically  using  a  series  of  short  laser  pulses  at  the  desired  sampling  rate.  These  pulses 
are  applied  in  parallel  to  the  modulators  a.'d  also  to  the  reference  channel. 
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Fig.  1-9.  (a)  Schematic  diagram  of  a  4-bit  electrooptic  A/D  converter, 

(b)  Intensity  vs  voltage  for  a  4-bit  A/D  converter  with  Gray-code  output 
(adapted  from  Ref.  25). 


*  1 

=  kL1V  , 

=  ir/2 

% 

=  2n"ZkL1V  , 

o 

II 

c 

n  =  2,  3. .  .N 

In  the  Gray-code  approach  shown  here,  the  modulator  electrode  lengths  differ  by  a  factor 
of  2  between  bit  channels  except  for  the  two  shortest  modulators  which  have  electrodes  of  equal 
length  L..  In  general,  the  electrode  length  Ln  for  the  nth  modulator  of  an  N-bit  converter  is 

Lr  =  2n"2L1>  where  n  =  2,  3. .  .N  and  n  =  2  corresponds  to  the  next-to-most  significant  bit  and 
n  =  N  corresponds  to  the  least  significant  bit  (LSB).  The  electrode  length  for  the  most  signifi¬ 
cant  bit  (MSB)  is  Lj.  The  net  phase  change  in  this  structure  for  the  various  channels  is 


(1-13) 


where  V  is  the  applied  voltage  and  k  is  a  constant  depending  on  the  wavelength,  the  waveguide 
spacing  and  the  crystal  electrooptic  parameters.  For  an  N-bit  device,  a  tt/2  phase  shift  in  the 
LSB  channel  corresponds  to  a  change  of  one  quantization  step  q  in  the  appL  d  voltage  and  re¬ 
sults  in  a  one-bit  change  of  the  output  code.  For  an  applied  signal  V&  =  VQ(1  +  sin  wt)  and 
2^  quantization  steps,  q  =  Vq/2^  1  and  kL^  =  ir/VQ.  The  static  phase  shift  =  ir/2  is  ob¬ 
tained  with  a  DC  bias.  The  Gray  code  minimizes  errors  because,  as  shown,  each  quantization 
step  only  results  in  a  change  in  one  bit  channel,  unlike  the  conventional  offset  binary  code  in 
which  multiple  bits  may  change  at  each  quantization  step.  In  addition,  the  Gray-code  output 
yields  one  more  bit  of  precision  than  the  offset  binary  code  for  a  given  set  of  electrode  lengths 
and  applied  voltage. 

The  electrooptic  technique  has  a  number  of  advantages  relative  to  other  conversion  methods 
First,  the  device  is  relatively  simple  since  only  one  modulator  and  detector  are  required  per 
bit.  An  important  characteristic  is  the  decoupling  of  the  analog  signal  from  the  optical  sam¬ 
pling  pulse.  This  should  eliminate  the  distortion  effects  common  to  conventional  devices  that 
utilize  diode -bridge  sampling  (or  tracking)  circuits.  The  electrooptic  technique  lends  itself  to 
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very  high  speed  operation  (>1  GHz)  with  very  short  sampling  pulses  ( ~  1 0  psec,  if  desirable) 
because  of  the  inherent  properties  of  mode-locked  lasers  which  can  be  used  to  generate  the 
sampling  pulse  train.  The  device  is  relatively  compact,  especially  if  one  of  the  recently  dem¬ 
onstrated  mode-locked  diode  lasers3*3,31  is  used  as  a  source.  Finally,  it  may  be  possible  to 
fabricate  a  monolithic  integrated  version  of  the  device  since  all  the  components  needed  have 
been  demonstrated  in  GaAs -based  semiconductors. 


2.  Design  Considerations 

An  analytical  study  was  performed  to  determine  the  fundamental  design  considerations  for 
this  electrooptic  device  as  a  high-speed  converter.  An  error  analysis  was  conducted  which 
examined  the  effects  of  optical  transit  time  through  the  modulator,  laser  pulse  width,  and  jitter. 
Analyses  were  also  conducted  to  determine  the  laser  optical  power  requirements  as  well  as  the 
requirements  on  the  analog  driver. 

Transit  time  limitations  were  computed  by  assuming  that  the  time -varying  analog  signal 
was  applied  simultaneously  to  the  entire  electrode  and  calculating  the  net  phase  change  (p  of  a 
wavefront  of  the  laser  pulse  as  it  traveled  the  length  of  the  electrode.  For  the  n1*1  modulator, 
characterized  by  a  transit  time  Tn,  and  an  analog  signal  V(t)  =  V  (1  +  sinwt), 

?  r  sin(i»)T  /2)  i 

^n-*2  [1+  (o>Tn/2)'  8inb»(t  +  Tn/2)|  •  U-14) 


It  can  be  seen  from  Eq.  (1-14)  that  the  net  phase  change  exhibits  a  nonuniform  attenuation 
and  a  delay  dispersion  between  bits.  Both  these  effects  are  caused  by  variation  in  electrode 
transit  times  between  bit  channels  and  lead  to  non-monotonicity  effects.  It  should  be  noted  that 
the  simple  averaging  of  the  light  pulse  over  the  electrodes  is  a  linear  operation  of  integration. 

It  can  be  completely  described  in  terms  of  a  linear  filter  and  thus  does  not  introduce  errors  in 
the  conversion. 

Non-monotonicity  due  to  nonuniform  attenuation  occurs  because  the  phase  attenuation  in¬ 
creases  the  required  voltage  to  produce  a  given  number  of  ir  phase  shifts.  It  can  be  shown  that 
the  effect  first  occurs  when  the  next-to-last  transition  for  the  LSB  (corresponding  to 

near-maximum  voltages)  occurs  at  a  larger  voltage  than  the  last  "0"-"l"  transition  in  the  next- 
to-LSB  [see  Fig.I-9(b)].  This  leads  to  the  requirement  that  T^  <  22  ^2/irfs,  where  fg  is  the 
maximum  signal  frequency  component  to  be  digitized. 

Delay  dispersion  is  a  potentially  much  more  severe  limitation.  By  a  similar  argument  to 
that  given  above,  non-monotonicity  occurs  for  T„  <  21  /irf  .  However,  this  effect  is  known 
and  can  be  compensated  for  by  introducing  an  appropriate  delay  between  bit  channels.  A  con¬ 
venient  way  to  do  this  is  to  align  all  the  modulator  electrodes  along  their  centerline  as  illus¬ 
trated  in  Fig.  I-9(a). 

Laser  pulse  width  effects  were  calculated  by  assuming  that  the  detector  output  was  deter¬ 
mined  by  the  average  intensity  over  the  entire  optical  pulse  width.  Transit  time  effects  were 
also  incorporated  in  the  analysis.  Assuming  Vg  =  V  (1  +  sinwt),  a  worst  case  of  the  high- 
frequency  zero-crossing  was  examined.  Both  a  square  and  a  Gaussian  pulse  shape  were  analyzed. 

For  the  ideal  case  of  a  square  wave  of  width  D,  one  obtains 


I 


n 


sin(k  co  D) 

— I — S-rf —  cos  fk  (w  T  -2)1 
k  u  D  n  s  n  1 

n  s 
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where  kn  =  kLnVQ/2.  A  nonlinear  error  occurs  when  the  sinx/x  terms  change  sign.  If  errors 

are  to  be  limited  to  the  LSB,  then  (fff  , 

s  -\  / pT 

A  similar  analysis  of  the  more  realistic  case  of  a  Gaussian  pulse  of  shape  ~e  '  leads  to 


In~cos[kn(u)sTn-2)]e 


,222 
-k  w  s 
n  s 


The  limit  on  laser  pulse  width  occurs  when  the  amplitude  is  too  small  for  detection.  This 

can  be  approximated  by  requiring  that  the  exponential  term  be  above  a  detection  limit  t] 

2  2  2 

{i.  e.,  exp[— kni»)  s  ]  >jj},  and  leads  to  a  full  width  half  maximum  s^^  for  the  Gaussian  pulse  of 


sl/2 


V (In  2)  (In  l/n  ) 

-  2.-R-3- 

ir  fg2 


The  limitation  in  the  jitter  of  the  onset  of  the  sampling  interval  is  set  by  consideration  of  the 
high-frequency  zero-crossing  of  the  analog  signal.  If  the  requirement  taken  is  that  the  voltage 
uncertainty  AV  be  less  than  quantization  step  size,  then  the  allowable  jitter  time  At  is  given  by 
At  =  (7rfg2N)‘1. 

The  necessary  laser  optical  power  was  estimated  by  defining  an  allowable  probability  of  a 
2-bit  error.  The  worst  case  occurs  if  the  LSB  is  very  near  a  threshold  crossing  and  the  next- 
to-MSB  (or  the  MSB)  is  also  near  a  threshold  crossing.  Sufficient  laser  power  should  be  used 
to  ensure  that  to  within  the  chosen  error  probability  the  more  significant  bit  is  detected  cor¬ 
rectly.  This  requires  that  igm  An  >  x(2^*  V,r).  where  igm  is  the  maximum  signal  photocurrent 
deviation  from  1^,  in  is  the  rms  noise  current,  and  x  is  set  by  the  desired  error  probability 
assuming  Gaussian  noise  statistics.  If  it  is  assumed  that  the  detectors  are  avalanche  photodiodes 
of  gain  G  and  responsivity  r  and  that  the  light  coupling  coefficient  into  the  waveguides  is  k. ,  one 
can  write  igm  =  P.k^rG,  where  Pj  is  the  required  incident  power.  This  leads  to  the  requirement 
that 

P.  >inx2N‘1Ak.rG  . 


The  analog  driver  requirements  were  computed  by  assuming  the  modulator  impedance  could 

be  modeled  as  a  capacitance  and  that  contact  resistance  was  negligible.  The  rms  signal  power 
2 

P  was  estimated  as  V  /2R,  where  R  is  the  source  impedance  and  2V  is  the  voltage  required 
a  N-l  °  o 

to  generate  2  ir  phase  shifts  in  the  LSB  modulator.  [Although  the  analog  signal  is  of  the  form 

V&  ~  VQ(1  +  sinwt),  the  DC  term  VQ  can  be  introduced  by  a  level  shift  and  therefore  does  not 

affect  the  driver  power  requirement.  ]  R  was  computed  assuming  that  the  RC  product  of  the 

modulators  and  source  would  correspond  to  a  3-dB  frequency  twice  f  .  That  is,  the  driver- 

modulator  would  be  no  worse  a  filter  than  other  components  in  the  signal  processing  system. 

Using  expressions  for  the  capacitance  and  electric  field  profile  of  coplanar  electrodes,32  it  is 

found  that 

Pa  =  ,rklk2  ln  (4W2  /W1  >  • 

Here  kj  and  k2  are  defined  by  the  coplanar  capacitance  and  the  electrooptic  coefficients, 
respectively,  X  is  the  optical  wavelength  and  Wj  and  W.,  are  the  spacings  between  the  inner  and 
outer  edges  of  the  coplanar  electrodes  along  the  waveguide.  To  minimize  the  necessary  power. 
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one  clearly  should  minimize  A  and  W^  and  maximize  Lj.  Another  requirement  on  the  analog 
driver  amplifier  is  that  its  linearity  not  deviate  by  a  quantization  step  (i.e.,  linear  to  one  part 
in  2N). 

Based  on  this  analysis  and  a  consideration  of  waveguide  materials  and  technology,  it  is  con¬ 
cluded  that  the  development  of  a  6 -bit  1-G S/sec  converter  should  be  feasible.  This  converter 
which  will  be  reported  on  in  the  future  would  require  an  array  of  LiTaO^  electrooptic  modulators, 
a  doubled  mode-locked  Nd:YAG  laser,  and  ~1  W  rms  of  analog  signal  power. 

F.  J.  Leonberger 

C.  E.  Woodward 

D.  L.  Spears 
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IL  QUANTUM  ELECTRONICS 


A.  REMOTE  SENSING  OF  CO  AND  NO  WITH  FREQUENCY-DOUBLED 

C02  LASER  RADIATION 

A  program  to  use  frequency-doubled  CC>2  laser  radiation  for  remote  sensing  of  CO  and  NO 
molecules  will  be  undertaken  in  the  new  roof  facility,  with  the  high -repetition  rate  mini-TEA 
laser  providing  the  C02  radiation  and  CdGeAs2  serving  as  the  frequency-doubling  crystal. 
Since  the  laser  is  line  tunable,  remote  sensing  is  limited  to  those  C02  laser  lines  whose  doubled 
frequency  coincides  with  an  absorption  line  of  the  molecules.  At  sea  level,  we  define  coincidence 
as  a  frequency  difference  An  <0.1  cm  *.  A  determination  has  therefore  been  made  of  all  such 
coincidences  between  1  to  0  absorption  bands  of  CO  and  NO  and  frequency -doubled  12C160  laser 
lines,  and  of  the  relative  merits  of  each  of  these  coincidences  for  remote  sensing  purposes.  In 
addition,  the  laser  pulse  energy  requirements  for  direct  detection  of  topographically  backseat - 
tered  radiation  at  the  preferred  frequencies  have  been  investigated. 

A  list  of  all  coincidences  falling  within  our  definition  is  given  in  Table  II-l.  To  determine 
which  of  these  coincident  lines  are  most  suitable  for  use  in  a  remote  sensing  system,  considera¬ 
tion  was  given  to  the  available  pulse  energy  of  the  mini-TEA  laser  at  each  line,  the  absorption 
cross  section  of  each  of  the  CO  and  NO  lines  and  the  atmospheric  transmittance  at  the  corre¬ 
sponding  frequencies.  The  NO  absorption  band  falls  within  a  spectral  region  of  strong  H20  ab¬ 
sorption  and,  as  a  result,  the  doubled  P(24)  line  of  the  10. 6- urn  band  is  the  only  frequency- 
doubled  laser  line  which  can  be  used  for  remote  sensing  of  NO  without  excessive  HzO  absorption. 
Either  the  P(24)  or  the  R(18)  lines  of  the  9.6-nm  band  can  be  used  for  monitoring  CO.  The 
frequency-doubled  P(24)  line  is  slightly  preferable,  and  further  considerations  are  limited  to 
that  line. 

The  energy  requirement  calculations  assume  topographic  backscatter  from  a  diffuse  reflec¬ 
tor  located  a  distance  R  from  the  receiver,  a  direct  detection  system  and  the  use  of  a  differen¬ 
tial  absorption  technique  in  which  the  reflected  response  at  the  molecular  absorption  (coincidence) 
frequency  is  compared  with  reflection  at  a  neighboring  nonabsorbing  frequency.  The  P(26)  C02 
laser  lines  of  the  10.6-  and  9.6-|im  bands  are  suitable  as  the  neighboring  nonabsorbing  lines  for 
NO  and  CO,  respectively. 

It  can  be  shown  that,  to  a  good  approximation,  the  difference  in  reflected  signals  6  from 
the  on-  and  off-coincidence  radiation  due  to  the  presence  of  an  absorbing  molecule  is  given  by 


where  K  is  the  receiver  efficiency,  p  the  reflectivity  of  the  topographic  reflector,  Pq(Ac)  the 
output  power  at  the  coincident  wavelength  Ac,  <rm(*c)  is  the  absorption  cross  section  at  A  ,  and 
n  is  the  pollutant  density  which  extends  over  a  distance  L.  T(A,  R)  is  the  two-way  atmospheric 
transmittance  over  a  range  R  at  wavelength  A  defined  by 

T(A,R)  =  exp[— 2/3m(A)R]  (H-2) 

where  0m  is  the  volume  extinction  coefficient  of  the  atmosphere  at  A.  The  values  of  a  and 
at  the  chosen  radiation  frequencies  have  been  determined  and  are  given  in  Table  II- 2  for  the 
standard  midlatitude  winter  and  midlatitude  summer  atmospheres.2 
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TABLE  II- 1 

SPECTRAL  COINCIDENCES  BETWEEN  NO  AND  CO  ABSORPTION  LINES 

AND  FREQUENCY-DOUBLED  C02  LASER  RADIATION 

Molecule 

Absorption 

Line 

Identification 

Absorption 

Line 

Frequency^ 
(cm  ’) 

COj  Laser 

Line 

Identification 

Doubled 

COj  Loser 

Frequency^ 

(cm  ') 

Av  | 

(2,,C02_,'mol) 

(cm  ’) 

NO 

R(73/2)3/5 

1977.2721 

+R(42) 

1977.29546 

0.0234 

R(71/2)3/2 

1975.2971 

+R(40) 

1975.24230 

-0.0548 

R(53/2)]/2 

1954.3952 

+R(22) 

1954.42922 

0.0340 

R(37/2)3/2 

1935.4951 

+R(8) 

1935.41584 

-0.0793 

R(19/2)|/2 

1909.1281 

+P(8) 

1909.09160 

-0.0365 

«’X2)1/2 

1881.0398 

+P(24) 

1881.09766 

0.0579 

CO 

R(6) 

2169.19834 

XR(30) 

2169.27032 

0.07198 

R(2) 

2154.59598 

X  R(1 8) 

2154.60498 

0.00900 

Q(11) 

2140.96109 

XR(8) 

2140.92456 

-0.03653 

P(14) 

2086.32231 

XP(24) 

2086.32650 

0.00419 

+  indicates  (10°0)  to  (00*1)  band;  X  indicates  (10*0)  to  (0280)  band. 

t  NO  frequencies  from  Laboratoire  de  Spectroscopie  Moleculaire,  Univ.  de  Paris  VI,  1973 
(unpublished).  CO  frequencies  evaluated  from  Dunham  coefficients  given  by  H.  Kildal,  R.S.  Eng, 
and  A.H.M.  Ross,  J.  Mol.  Spectrosc.  53,  479  (1974). 

t  COo  laser  frequencies  from  Handbook  of  Lasers,  R.  J.  Pressley,  Ed.  (Chemical  Rubber  Company, 
Cleveland,  Ohio,  1971). 
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TABLE  11-2 

FREQUENCY-DOUBLED  LINES  FOR  REMOTE  SENSING 

COj  Loser 

Doubled 

Radiated 

Absorption 
Cross  Section 

o  (cm  ) 
m 

Volume  Extinction 

Coefficient  (3  (km  ' ) 
m 

Line 

Identification 

Frequency 
(cm  ') 

Mid  latitude 
Winter 

Midlatitude 

Summer 

P<24) 

(10°0)-  (02°0) 

2066.3265 

9.08X  KT19 

0.1654 

0.292 

P(26) 

(1 0°0)  —  (Q2°0) 

2082.5582 

0.0334 

0.1077 

P(24) 

(10°0)-  (00»1) 

1881.0977 

4.85  X  10-20 

0.335 

1.445 

P(26) 

(10°0)-  (00°1) 

1877.3780 

0.411 

1.694 

The  minimum  detectable  value  of  6  is  defined  to  occur  for  a  signal -to -noise  ratio  S/N  =  1; 
i.e.,  for  6  =  PN  =  (NEP)  (2Af)*  ,  where  the  noise  is  taken  as  dark-current  limited  in  the  in¬ 

frared.  The  minimum  energy  per  pulse  is  then 


.min 
o  ’  c' 


p  min 

o 


(*c>  T 


2irR2  (NEP)  s/T _ 

KpA  T(Xc)  {1  -  exp[-2nam(Xc)  L)} 


(II- 3) 


where  we  have  taken  the  time-bandwidth  product  rAf  =  2.  When  averaged  over  a  train  of  N  pul¬ 
ses,  the  minimum  energy  requirement  is  reduced  by  the  factor  (N)  To  evaluate  Eq.  (II- 3), 

the  values  NEP  =  3  x  10  Hz  p  =  0.5  and  K  =  0.1  are  assumed?  Our  effective  telescope 

2  -7 

area  A  “  600  cm  and  t  =*  10  sec.  Then 


E  (A  )  >. 
o  c'  x 


_ 2  k  10~6  R2 _ 

T(*c)  •m  fl-exp(-2n<7m(Xc)L)} 


(II- 4) 


The  minimum  energy  required  per  pulse  is  evaluated  as  a  function  of  range  in  Fig.  II- 1  for 
a  CO  concentration  of  100  ppm  over  a  distance  L  =  50  m  for  midlatitude  winter  and  summer  at¬ 
mospheres.  The  minimum  energy  vs  range  curves  for  a  similar  concentration  of  NO  are  given 
in  Fig.  II- 2. 

Using  a  mini-TEA  C02  laser  with  a  "three-mirror"  grating-tuned  cavity,4  second -harmonic 
outputs  of  over  2  mj/pulse  have  been  obtained  from  the  P(24)  lines  of  both  the  9.6-  and  10.6-um 
bands.  A  repetition  rate  of  500  Hz  has  been  achieved  with  the  mini-TEA  laser;  thus  averaging 
pulses  over  a  one-second  interval  can  lead  to  an  effective  pulse  energy  in  excess  of  40  mj.  Lon¬ 
ger  averaging  times  would  further  increase  the  effective  range  or  S/N  ratio. 

The  above  results  indicate  that  remote  sensing  of  CO  and  NO  from  topographic  targets  over 
ranges  of  a  few  kilometers  should  be  feasible  with  our  system.  The  range  is  most  limited  for 
NO  under  conditions  of  high  atmospheric  water  vapor  content  due  to  the  large  H20  absorption  in 
the  NO  absorption  band  region.  Although  the  standard  midlatitude  atmospheres  encompass  the 
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Fig.  II- 1.  Minimum  energy  per  pulse 
required  for  S/N  =  1  as  a  function  of 
range  assuming  a  CO  concentration 
of  100  ppm  over  a  distance  L  =  50  m 
for  standard  midlatitude  summer  and 
winter  atmospheres  and  topographic 
reflection. 


Fig.  II- 2.  Minimum  energy  per  pulse 
as  a  function  of  range  assuming  an  NO 
concentration  of  100  ppm  over  a  dis¬ 
tance  L  =  50  m  for  standard  midlati¬ 
tude  summer  and  winter  atmospheres 
and  topographic  reflection. 
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water  vapor  normally  encountered,  an  accurate  determination  at  the  time  of  the  experiment  will 
be  required.  In  addition,  the  level  of  ambient  CO  must  be  established  for  CO  monitoring.  The 


standard  atmospheres  assume  a  CO  content  of  0.07  5  ppm  (Ref.  2).  This  is  somewhat  lower  than 
the  global  average,5  and  is  much  lower  than  values  which  have  been  measured  in  urban  environ¬ 
ments.5  The  presence  of  additional  amounts  of  ambient  CO  strongly  increases  the  opacity  of  the 
atmosphere  at  the  coincident  frequency.  For  example,  an  increase  in  ambient  CO  content  to 
0.3  ppm  will  result  in  a  tenfold  decrease  in  atmospheric  transmittance  over  a  2-km  range,  with 
a  corresponding  increase  in  the  minimum  energy  required.  Finally,  the  broad  gain- bandwidth 
of  our  frequency- doubled  atmospheric  laser  will  alter  the  effective  value  of  am  from  the  values 
shown  in  Table  II- 2,  which  are  based  on  the  line  intensities  given  by  McClatcheyet  al.2  and  the 
frequency  difference  Av.  To  establish  the  proper  value  for  our  measurements,  it  will  therefore 
be  necessary  to  determine  the  effective  values  of  am  for  CO  and  NO  using  the  frequency-doubled 


mini-TEA  laser  as  the  coherent  radiation  source. 


N.  Menyuk 


B.  THERMAL  BEHAVIOR  OF  NdP50J4  LASERS 

During  experimental  investigations  of  high-Nd-concentration  lasers,  particularly 
NdP,0,  .(NPP),  it  has  been  observed  that  the  output  power  often  decreases  at  high  excitation 

D  14  7  g 

levels  or  duty  cycles.  '  In  order  to  develop  useful  CW  NPP  lasers  it  is  necessary  to  understand 
the  mechanisms  which  cause  this  effect. 

A  list  of  possible  mechanisms  for  thermally  caused  degradation  of  the  laser  performance 
includes: 

4 

(1)  increased  population  of  the  Ij^/2  l°wer  laser  level,  reducing  the 
population  inversion  and  gain; 

4 

(2)  redistribution  of  population  in  the  F^2  levels,  reducing  the  gain  on  the 
dominant  Rj-Yj  laser  transition; 

(3)  increase  in  laser  transition  linewidths; 

(4)  decrease  in  fluorescence  lifetime; 

(5)  thermal  lensing,  through  direct  temperature  dependence  of  the  index  of 
refraction  or  optical  path  length,  or  indirect  strain  and  elasto-optic 
effects; 

(6)  optical  deformation  from  ferroelastic  domains. 

The  main  conclusion  we  have  reached  is  that  for  our  NPP  laser,  direct  lower  laser  level  popula- 
34  4 

tion  in  the  Nd  ln/2  level  (item  1)  is  the  dominant  cause  of  efficiency  decrease. 

A  schematic  drawing  of  the  laser  apparatus  is  shown  in  Fig.  II-  3.  The  NPP  laser  had  mir¬ 
rors  with  5-cm  radii  of  curvature.  One  mirror  had  high  reflectivity  at  1052  nm  and  the  other 
had  reflectivity  of  either  99.6  or  99.0  percent.  The  laser  rod  in  the  center  of  the  cavity  was 
2  mm  square  and  5  mm  long  along  the  lasing  ('c')  axis.  The  end  faces  were  polished  parallel  to 
20  arc-sec  but  were  not  antireflection  (AR)  coated.  The  effective  mirror  separation  for  mode 
calculations  is  Leff  =  LM  —  f  c(n  -  l/n),  where  LM  is  the  actual  mirror  separation,  lc  is  the 
laser  rod  length  (5  mm),  and  n  is  its  index  of  refraction  (1.60).  In  our  cavity  Le^  =  9.5  cm, 
which  is  sufficiently  far  removed  from  the  concentric  case  (L^  =  2  x  mirror  radius  =  10  cm) 
that  the  cavity  should  be  quite  stable  to  perturbations. 
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The  NPP  laser  was  pumped  with  collinear  excitation  from  a  CW  rhodamine  6G  dye  laser 
tuned  to  590  nm.  This  wavelength  was  chosen  by  tuning  the  dye  laser  until  the  pump  beam  was 
barely  transmitted  by  the  crystal,  so  as  to  obtain  nearly  complete  absorption  as  uniformly  as 
possible  within  the  crystal.  A  7.5-mm-dia.  aperture  was  put  in  the  beam  before  the  pump  focus¬ 
ing  lens  and  Si  photodiode  used  to  measure  the  pump  power.  This  diode  was  calibrated  to  mea¬ 
sure  the  590  nm  power  transmitted  through  the  end  mirror  (high- reflectance  at  10  52  nm).  The 
laser  output  power  was  measured  with  a  calibrated  Si  photodiode  preceded  by  a  filter  to  eliminate 
the  slight  amount  of  transmitted  pump  radiation.  A  Gian  prism  could  be  inserted  before  the  pho¬ 
todiode  to  measure  the  laser  polarization.  To  examine  the  laser  spectrum,  the  laser  output 
beam  was  transmitted  by  a  fiber  optic  bundle  to  a  1-m  double-grating  spectrometer,  followed 
by  a  cooled  InGaAs  photomultiplier,  boxcar  amplifier,  and  chart  recorder.  The  laser  mode 
profile  was  obtained  with  an  infrared -sensitive  Si  TV  vidicon.  Thus,  in  our  series  of  experi¬ 
ments  we  measured  the  input  power  and  duty  cycle,  and  the  output  power,  spectrum,  polariza¬ 
tion,  and  mode  profile.  In  the  data  below,  the  measured  powers  refer  to  the  quasi-CW  powers 
during  the  excitation  pulse,  not  the  average  powers. 

Figure  II- 4  shows  the  experimental  data  for  the  output  vs  input  power  with  the  99.6-percent- 
reflecting  (0.4-percent-transmitting)  mirror  using  duty  cycles  of  rj  =  0.25,  0.50,  and  1.0.  With 
increasing  duty  cycle  and  heating,  the  threshold  increases  and  the  output  and  efficiency  decrease. 
During  the  course  of  these  measurements  the  laser  mode  profile  was  measured,  and  no  change 
in  shape  or  mode  waist  was  observed.  Data  similar  to  Fig.  II- 4  were  obtained  with  the 
99.0-percent-reflecting  mirror  and  are  shown  in  Fig.  II- 5.  The  larger  output  transmission  causes 
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Fig.  II- 4.  Output  power  as  a  function 
of  input  power  for  NPP  laser  with 
0.4-percent-transmitting  output  mir¬ 
ror.  The  different  data  and  curves 
refer  to  duty  cycles  of  0.25,  0.50, 
and  1.0. 


Fig.  II-  5.  Output  power  as  a  function 
of  input  power  for  NPP  laser  with 
1.0-percent-transmitting  output  mir¬ 
ror.  The  different  data  and  curves 
refer  to  duty  cycles  of  0.25,  0.50, 
and  1.0. 
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somewhat  higher  thresholds,  but  also  gives  greater  slope  efficiencies.  Heating  effects  are  again 
evident.  From  these  data,  the  internal  loss  is  calculated  to  be  0.67  percent. 

After  performing  these  measurements,  the  laser  rod  was  remounted  for  better  cooling  with¬ 
out  changing  the  mirror  separation.  The  'a'  side  faces  of  the  rod  were  held  between  two  small 
thermoelectric  cooler  modules.  A  lowering  of  the  surface  cooling  temperature  of  about  1 2*C 
restored  the  high -duty -cycle  performance  to  CW  operation. 

On  the  basis  of  calculations  and  our  experimental  results,  we  have  found  that  the  terminal 
state  thermal  population  is  the  most  important  effect  of  those  listed  above;  to  verify  this  conclu¬ 
sion  further,  we  have  calculated  its  contribution.  The  terminal  state  population  has  an  exponen¬ 
tial  temperature  dependence.  Using  measured  values  of  the  thermal  conductivity  in  NPP,  a 
simple  model  was  used  to  calculate  the  temperature  rise  in  the  laser  mode  volume,  leading  to 
a  nonlinear  equation  for  the  resonant-loss  threshold  term,  proportional  to  the  terminal  state 
population.  The  results  of  this  calculation  are  shown  as  the  curves  in  Figs.  II- 4  and  -5.  A  heat 
sink  impedance  parameter  was  adjusted  in  order  to  fit  the  data.  In  view  of  the  approximations 
in  the  model,  the  experimental  agreement  is  quite  reasonable,  and  shows  that  this  simple  in¬ 
trinsic  heating  effect  can  explain  the  observed  power  dependence  of  the  laser  output. 

While  these  results  pertain  to  one  particular  example,  the  general  method  and  conclusions 
should  have  a  considerably  wider  range  of  applications.  For  example,  a  useful  geometry  for 
practical  NPP  laser  devices  is  the  transverse  pumping  arrangement,  using  either  CW  sources 
such  as  light-emitting  diodes  (LEDs)9  or  pulsed  sources  such  as  flash  lamps.10  In  such  cases 
the  pump  absorption  will  be  even  more  uniform  than  in  the  focused  collinear  pumping  method, 
and  our  conclusions  should  be  valid  for  the  transverse  pumping  schemes  as  well.  Generally,  in 
such  devices  the  geometric  pumping  efficiency  will  be  lower  than  in  the  collinear  case,  and  a 
larger  fraction  of  the  pump  will  heat  portions  of  the  laser  rod  outside  the  laser  mode  volume. 
This  should  affect  only  the  calculation  of  the  linear  heating  coefficients.  The  only  restriction  to 
our  conclusions  that  increased  population  of  the  levels  is  the  most  significant  thermal 

effect  in  NPP  lasers  is  that  the  mode  be  stable  with  respect  to  small  thermo-optical  effects. 

If  heating  is  a  problem  in  CW  or  high -repetition- rate  NPP  lasers,  different  remedies  are 
possible.  First,  active  cooling  of  the  laser  rod  may  be  sufficient  to  lower  the  temperature  in 
the  lasing  volume.  Second,  a  lower  Nd  concentration  would  decrease  both  the  population 

and  the  nonradiative  decay  rate,  although  the  pump  absorption  coefficient  would  also  decrease. 
Large,  high-quality  laser  rods  of  Nc^Laj _XP501 4  (x  <*  0.25  to  0.75)  have  been  difficult  to  fabri¬ 
cate  because  of  crystal  strain,  but  this  problem  eventually  may  be  overcome,  or  other  high-Nd- 
concentration  materials  might  be  used.  Finally,  it  may  be  desirable  to  use  the  1320-nm  transi¬ 
tions  from  F3/2  ^13/2-  wh°se  terminal  state  population  is  negligible.  Although  the  gain 

cross  section  is  ~5  times  smaller  than  at  10  50  nm  (Refs.  11  and  12),  resonant  absorption  and  its 

thermal  dependence  should  not  be  important.  „  _  . 

S.  R.  Chinn 
W.K.  Zwicker^ 

C.  SPECTROSCOPIC  MEASUREMENTS  OF  Tb^j _xP5Oj4  COMPOUNDS 

A  series  of  measurements  has  been  initiated  to  determine  the  potential  of  terbium -yttrium 
pentaphosphate  ((Tb,  Y)PP]  crystals  as  laser  materials  for  the  blue-green  spectral  region.  Ab¬ 
sorption  and  emission  spectra  have  been  taken  and  key  radiative  and  nonradiative  rates  have 


t  Philips  Laboratories,  Briarcliff  Manor,  New  York. 
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been  measured  with  pulsed  rare-gas:halide  laser  excitation.  A  strong  absorption  band  in  coin¬ 
cidence  with  the  KrF  pump  laser  has  been  found  corresponding  to  the  lowest  4f-5d  transition. 

This  may  provide  an  excellent  way  in  which  to  pump  the  candidate  4f  transitions  since  5d  —  4f 
transfer  is  observed  to  be  rapid. 

5  5 

In  100%  TbPP  the  decay  times  of  the  prospective  upper  laser  levels  (  D3  and  D4)  are 
4.1  nsec  and  3.1  msec,  respectively,  with  the  5D,  decay  occurring  predominantly  into  the  5D. 

^  5  7  4 

state  by  nonradiative  processes.  The  strongest  4f  —  4f  radiative  transitions  are  D.  —  Fc  and 

C  7  4  3 

"X)3  F^.  Growth  techniques  developed  for  NdPP  crystals  appear  to  be  applicable  to  the  growth 

of  (Tb,  Y)PP  materials  with,  however,  some  alteration  in  growth  parameters  due  to  the  difference 
between  Tb,  and  Nd  ionic  radii.  Based  on  the  measured  decay  times  and  transition  bandwidths, 
it  is  estimated  that  oscillators  can  be  constructed  from  these  materials  with  rare-gas:halide 
laser  pumping  if  adequate  crystals  of  high  Tb  doping  can  be  obtained. 

The  (Tb,  Y)PP  materials  are  also  good  systems  in  which  to  study  both  5d  —  4f  and  4f  — 4f 

5  5 

energy  transfer.  A  particularly  interesting  example  is  the  D3  —  D4  transfer,  which  presum¬ 
ably  occurs  as  a  result  of  long  range  interionic  forces. 

D.  J.  Ehrlich  R.  M.  Osgood,  Jr. 

S.  R.  Chinn  W.  K.  Z  wicker  t 


D.  CW,  HIGH- POWER  Ni:MgO  LASER 

The  first  CW  operation  of  a  Ni:MgO  solid-state  laser  has  been  observed.  Nearly  6  W  of  out- 

3  3 

put  at  1.32  um  was  obtained  by  optical  pumping  of  the  A ^  —  T  ^  transition  in  a  1.85-cm-long 
crystal  of  Ni:MgO  cooled  in  a  liquid  nitrogen  Dewar.  Up  to  30  percent  of  the  absorbed  1.06-|im 
pump  radiation  from  a  Nd:YAG  laser  was  converted  to  Ni:MgO  laser  output.  Quasi-CW  opera¬ 
tion  was  observed  up  to  a  temperature  of  235  K.  Preliminary  temperature  tuning  results  showed 
operation  in  several  wavelength  regions  between  1.316  and  1.409  um. 

The  laser  crystal  was  grown  by  a  carbon-arc  fusion  technique  developed  and  carried  out  at 

1 3 

Oak  Ridge  National  Laboratory.  5  wt%  of  NiO  was  added  to  the  starting  growth  material;  the 
resultant  crystals  contained  about  0.5  wt%  of  Ni.  The  crystal  used  in  the  experiments  was 
1.85  cm  long  and  had  a  cross  section  of  0.4  x  0.5  cm.  The  ends  of  the  crystal  were  polished  flat 
and  AR  coated  for  1.4  um. 

Laser  operation  was  observed  with  the  crystal  in  the  center  of  a  2  5-cm-long  confocal  cavity 
with  1.2-percent  output  coupling  at  1.32  um.  The  crystal  was  clamped  to  a  heat  sink  attached  to 
the  cold  finger  of  a  liquid  nitrogen  Dewar.  The  heat  sink  temperature  could  be  raised  above  77  K 
by  heating  elements  mounted  in  the  heat  sink  and  monitored  by  a  Pt  resistor.  The  two  laser  mir¬ 
ror  surfaces  were  inside  the  vacuum  chamber  of  the  Dewar.  The  pump  laser  output  was  focused 
in  the  crystal  by  a  15-cm  focal  length  lens.  About  50  percent  of  the  incident  radiation  was 
absorbed. 

A  curve  of  Ni:MgO  laser  output  vs  absorbed  pump  power  input  is  shown  in  Fig.  II- 6.  The 
heat  sink  temperature  varied  from  around  82  K  at  the  lowest  pump  power  to  95  K  at  22.6  W  of 
absorbed  power.  At  95  K  the  temperature  was  fluctuating  rapidly,  presumably  because  of  rapid 
boiling  at  the  liquid -nitrogen  cold-finger  interface.  Higher  output  powers  should  be  obtainable 
with  an  improved  cold-finger  design.  The  pump  laser  was  highly  multimode  (>10  times  the  dif¬ 
fraction  limit)  while  the  output  of  the  Ni;MgO  laser  for  the  curve  shown  in  Fig.  II- 6  was  a  mixture 

t  Philips  Laboratories,  Briarcliff  Manor,  New  York. 
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Fig.  II- 6.  CW  output  power  of  Ni:MgO 
laser  vs  absorbed  1.06-nm  pump  power. 
Maximum  conversion  efficiency  of  pump 
into  laser  power  ~30  percent.  Roll-off 
in  output  power  was  due  to  limited  cool¬ 
ing  capacity  of  LN2  Dewar. 


Po».  <"> 

of  the  TEM00  and  TEMQ1  modes.  Such  behavior  is  possible  because  the  cavity  TEMQ0  mode 
size  in  the  Ni:MgO  crystal,  230  11m,  is  ~5  times  larger  than  the  expected  size  of  a  focused 
TEMqq  mode  from  the  pump  laser.  By  reducing  the  number  of  modes  operating  in  the  pump 
laser  it  was  possible  to  obtain  single  TEMQq  mode  output  from  the  Ni:MgO  laser,  with  2. 5  W  of 
CW  power. 

Preliminary  results  of  the  effects  of  heat  sink  temperature  on  laser  performance  indicated 
that  laser  operation  was  possible  up  to  235  K  at  a  50-percent  duty  cycle  with  a  6-times-threshold 
increase  over  the  threshold  at  82  K.  As  the  temperature  increased  from  82  K,  the  operating 
wavelength  tuned  from  1.316  to  1.328  nm  at  131  K,  shifted  to  1.369  nm  at  153  K  and  further  shifted 
to  1.409  11m  at  235  K. 

More  detailed  measurements  on  the  temperature  behavior  of  the  Ni:MgO  laser  will  be  per- 

14 

formed.  Since  the  lifetime  of  the  laser  transition  remains  constant  up  to  300  K,  it  is  possible 
that  room  temperature  operation  of  the  Ni:MgO  laser  can  be  obtained. 

P.  F.  Moulton 
A.  Moo  radian 
Y.  Chent 

E.  MEASUREMENT  OF  THE  ELECTRONIC  THIRD-ORDER  SUSCEPTIBILITY 
IN  LIQUID  CO-Oz 

We  have  measured  the  electronic  contribution  to  the  nonresonant  third-order  suscepti- 
( 3)  ei 

bility  x  of  a  liquid  CO-O,  mxture.  This  result  is  important  for  determining  the  optimum  con- 
nr  *  1516 

ditions  for  the  two-photon  resonant  third-harmonic  generation  (THG)  process  ’  in  liquid 

CO-O,  and  is  also  of  fundamental  interest. 

C  (3)  (3) 

In  addition  to  x  \  ,  there  are  also  contributions  to  x  due  to  electrostriction  and  molecular 
el  ,  nr 

reorientation.  These  are  larger  than  at  low  driving  frequencies,  but  make  negligible  con¬ 
tributions  to  frequency  mixing  processes  at  optical  frequencies  because  of  their  limited  frequency 

( 3) 

response.  A  THG  technique  is  therefore  well  suited  for  measuring  X^j  •  Direct  measurements 
t  Oak  Ridge  National  Laboratory. 
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are  difficult  because  of  calibration  problems;  we  have  therefore  chosen  an  interference  technique. 

By  tuning  a  CO,  laser  through  the  vibrational  two-photon  resonance  in  liquid  CO-O,,  a  variation 

£  ( 3)  £ 

in  the  THG  signal  is  observed  due  to  the  interference  between  \  ,  '  and  the  two-photon  resonant 

( 3)  el 

susceptibility  x  By  measuring  the  frequency  difference  between  maximum  and  minimum 

res  (3)  (3) 

third -harmonic  signal,  \  \  can  be  determined  relative  to  x  '• 

, ,,  ei  res  17-21 

The  total  x  '  for  liquid  C0-02  mixtures  that  contributes  to  THG  is  of  the  form 


x(3)  (a;,  Ci),  w)  =  ^  Aw  I  x  ^ I  i 
*  w,  w,  w,  2  “"nIXreslpk 


(II- 5) 


in  the  limit  of  wtq  »  1,  where  w  is  the  width  at  77  K  of  the  O  and  S  bands  for  gaseous  CO 
vibrational  two-photon  absorption  and  is  the  CO  dielectric  relaxation  time  in  liquid  CO-O,. 

V7  t 

There  are  two  contributions  to  the  two-photon  resonant  susceptibility  with  very  different  line- 
widths  Awn  and  Aw^.  The  narrow  (Awn)  and  broad  (Aw^)  resonance  components  are  associated 
with,  respectively,  the  isotropic  or  scalar  part  a  and  the  anisotropic  part  p  of  the  CO  vibra¬ 
tional  Raman  polarizability  tensor.1^’20  In  Eq.  (II- 5)  it  is  assumed  that  4y  2  Awn/45  Aw^  «  1, 
where  y  =  p/a  =  1.53  for  CO  (Ref.  22).  The  linewidth  and  the  vibrational  frequency  wQ  of  the  nar¬ 
row  resonance  component,  including  the  CO  concentration  dependence,  have  been  measured  by 
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polarized  spontaneous  Raman  scattering.  The  broad  linewidth,  given  by  Aw,  *=  6/r_  K  60  cm 

'  17  24  d  tf 

can  be  estimated  from  infrared  absorption  data.  ’  The  center  frequency  of  the  broad  reso¬ 
nance  differs  slightly  from  the  position  of  the  narrow  component.  The  difference  is  negligible, 
however,  compared  to  Aw^. 

Figure  II- 7  shows  the  frequency  dependence  of  the  observed  third-harmonic  signal  for  a 

CO-O,  mixture  at  77  K  with  a  0.077  relative  CO  peak  absorbance.  This  relative  peak  absorbance 
£  17 

is  approximately  proportional  to  the  relative  CO  concentration.  The  measurements  were  car¬ 
ried  out  using  a  tight- focusing  geometry  with  a  0.65-cm  confocal  parameter  in  a  9.6-cm  cell. 

The  changes  in  the  phasematching  condition  and  the  confocal  parameter  were  negligible  for  the 
frequency  range  of  1035.46  to  1084.62  cm  1  over  which  the  CO,  laser  was  tuned  stepwi  The 
third -harmonic  signal  in  Fig.  II-7  is  therefore  directly  proportional  to  |  x'  (oi,  w,  w)  |  .  The  sig¬ 
nal  varies  by  more  than  six  orders  of  magnitude.  A  destructive  interference  between  the  reso- 
nant  and  nonresonant  terms  in  x  (w,  w,  w)  is  evident  on  the  high-frequency  side  of  the  two-photon 
resonance  at  a  frequency  mismatch  6w  =  wQ  —  2w  of  about  —6  cm  *.  Minimum  signal  is  obtained 


Fig.  II-7.  Interference  between  resonant  and 
nonresonant  contributions  to  x  ^  3Kw,  w,  w)  for 
aCO-02  mixture  with  0.077  CO  relative  peak 
absorbance. 
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( 3) 

when  the  real  part  of  x  (w,  w,  o>)  vanishes.  This  gives  a  relationship  between  the  nonresonant 
and  resonant  parts  of  the  third-order  susceptibility  tensor  given  by  «  —  Aw>n/ 26w  Ix^es lpk 

when  a  small  correction  due  to  the  broad  resonance  is  neglected.  The  third-harmonic  signal 
does  not  vanish  at  maximum  destructive  interference.  The  residual  signal  is  determined  mainly 
by  the  imaginary  component  of  the  broad  resonance. 

In  order  to  determine  accurately  the  signal  level  and  frequency  position  for  maximum  de¬ 
structive  interference  we  have  tuned  through  this  region  by  changing  the  CO  concentration  while 
keeping  the  pump  frequency  constant  at  1071.88  cm  *  [CO^  R(10)  line].  Figure  II- 8  shows  the 
measured  concentration  dependence  of  the  third -harmonic  signal.  A  minimum  is  observed  at  a 
CO  relative  peak  absorbance  of  0.081.  Without  the  interference,  the  third -harmonic  signal  would 
have  exhibited  a  single  peak  characteristic  of  a  tight -focusing  geometry. 


Fig.  II- 8.  Third -harmonic  signal  near 
cancellation  of  resonant  and  nonresonant 
contributions  to  x^3Hw,  to,  w)  in  liquid 
CO-O2  (laser  emission  at  1071.88  cm-1, 
confocal  parameter  is  0.65  cm). 


The  solid  curves  in  Fig.  II- 7  and  -8  are  the  results  of  a  computer  calculation  using  x 


el3>'  A“b* 


and  the  signal  amplitudes  as  the  only  adjustable  parameters.  In  Fig.  II- 8  the  dependence  of 
(3) 

X  .  and  Aw.  on  the  CO  concentration  was  neglected,  and  the  phasematching  integral  was  calcu- 

^  25  (3) 

lated  using  the  wave  vector  mismatch  reported  previously.  The  value  of  Xgj  is  most  sensitive 

to  the  curve  fitting  in  Fig.  II -  8.  We  obtain 


'  (3>  I  u/x  V*  =  63  ±  2 

'res  pk  Ael 


(3) 


for  a  0.077  CO  relative  peak  absorbance.  Using  our  previously  reported  value  5  for  |x  |  .  , 

(3)  -34  /  3  res  pK  -1 

we  obtain  Xgj  =  3.9  x  10  Asm/V  .  The  broad  resonance  has  a  fitted  linewidth  of  Au^  =  62  cm 

which  is  consistent  with  the  estimated  value. 

It  is  interesting  to  examine  the  validity  of  applying  Miller's  rule2^  to  the  third -order  elec¬ 
tronic  susceptibility  in  liquids.  The  electronic  susceptibility  for  CS,  has  been  reported27  as 


=  2. 2  x  10  32  Asm/V3.  Using  the  indices  of  refraction  n_„  =  1.62  and  n 

Lb,  ' 


xel,CS, 


we  obtain  from  Miller's  rule 


CO-O. 


=  1.23, 


nco-o2  1 


“cs. 


- 1 


4 

^  xel,CS2  =  2’2  X  10'34  Asm/v3 
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which  is  in  reasonable  agreement  with  the  measured  value.  The  indices  of  refraction  are  about 

the  same  for  CO  and  O,  liquids.  Therefore,  from  this  result,  we  expect  little  variation  in  mea- 
1 3)  * 

sured  x  with  CO  concentration.  H  Kildal 

S.  R.  J.  Brueck 

F.  EFFICIENT  INFRARED  AC  KERR  SWITCHES  USING  SIMPLE 
CRYOGENIC  LIQUIDS 

We  have  constructed  efficient  infrared  AC  Kerr  switches  of  visible  radiation  using  the  cryo¬ 
genic  liquids  O^,  N^,  and  CO  as  the  Kerr-active  media.  An  intense  CO^  TEA  laser  pulse  has 
been  used  to  induce  a  birefringence  greater  than  i/2  wave  for  a  visible  probe  beam  in  liquid  O^. 
Over  67  percent  of  the  visible  laser  power  was  switched  into  the  orthogonal  linear  polarization 
for  a  CO^  input  power  of  1  MW.  These  infrared -visible  Kerr  switches  offer  significant  promise 
for  the  generation  of  subpicosecond  infrared  pulses  using  an  intense  visible  laser  pulse  to  open 
the  switch  and,  conversely,  for  upconverting  the  temporal  dependence  of  short  C02  pulses  into 
the  visible  spectral  region  with  subpicosecond  resolution  for  diagnostic  purposes. 

These  cryogenic  liquids  are  transparent  throughout  the  infrared  with  the  exception  of  a  few 
relatively  narrow  absorption  bands  corresponding  to  the  vibrational  fundamentals  and  overtones 
of  the  diatomic  molecules.  The  dominant  contribution  to  the  Kerr  coefficient  arises  from  the 
partial  orientation  of  the  anisotropic  molecules  in  the  intense  field;  the  speed  of  response  of  the 

medium  is  set  by  the  orientational  relaxation  time  r  which  is  in  the  0.1-  to  0.3-psec  range  for 
24  28  ^ 

the  cryogenic  liquids.  ’  This  is  substantially  shorter  than  the  relaxation  time  for  the  more 
conventional  Kerr  liquids  CS.,  (rr  =  1.8  psec)  (Ref.  29)  and  nitrobenzene  (rr  =  32  psec)  (Ref.  30) 
which  have  been  used  in  the  visible  spectral  region.  Due  to  their  infrared  absorption,  these 
more  complex  liquids  are  not  generally  useful  throughout  the  infrared. 

Efficient  Kerr  switches  can  be  used  to  upconvert  the  temporal  information  on  an  infrared 
laser  pulse  to  the  visible  spectral  region.  Conversely,  these  Kerr  switches  along  with  pulsed 
visible  and  near- infrared  lasers  may  be  used  to  gate  infrared  laser  pulses.  The  required  laser 
powers  are  in  the  1-  to  10-MW  range.  In  order  to  realize  the  material  limited  switching  speeds 
of  0.1  to  0.3  psec  it  will  be  necessary  to  compensate  for  the  difference  in  the  pulse  group  veloc¬ 
ities  between  the  infrared  and  optical  spectral  regions.  This  can  be  achieved  by  adding  a  mole¬ 
cule  which  exhibits  the  appropriate  dispersion  to  the  liquid,  as  has  recently  been  demonstrated 
for  phasematching  of  the  third -harmonic  generation  process  in  liquid  CO-O^  mixtures.  ^  The 
dispersion  of  all  these  liquids  is  sufficiently  small  that  pulse  spreading  of  the  individual  pulses 
is  negligible  for  0.1-psec  duration  pulses  and  10-cm  path  lengths. 

The  details  of  the  Kerr  switches  and  the  measurement  of  the  AC  Kerr  coefficients  of  the 
cryogenic  liquids  will  be  reported  elsewhere. 

S.  R.  J.  Brueck 
H.  Kildal 

G.  LASER  PHOTODEPOSITION  OF  METAL  FILMS  WITH  MICROMETER 
SIZE  FEATURES 

We  have  demonstrated  the  use  of  both  pulsed  and  CW  UV  lasers  to  deposit  metal  structures 
with  micrometer  size  features.  Our  experiments  have  explored  one  example  of  a  general  class 
of  laser-initiated  heterogeneous  photochemical  reactions  which  have  application  to  removal  of 
material  (etching)  as  well  as  to  the  deposition  of  elements  or  compounds.  The  use  of  the  laser 
permits  ready  spatial  localization  of  the  reaction  region;  hence,  laser  photodeposition  may  be 
applicable  to  the  production  of  structures  useful  for  microelectronics. 
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These  techniques  should  be  distinguished  from  laser  processing  based  on  strictly  thermal 
effects  such  as  laser  annealing  and  chemical  vapor  deposition  using  laser  heated  surfaces.  In 
our  experiments  molecular  bonds  are  broken  directly  by  UV  laser  photons,  and  free  metal  atoms 
are  produced.  These  atoms  then  condense  on  a  substrate. 

We  used  the  metal  alkyl  compounds  trimethylaluminum  (TMA1),  A1(CH3>3,  and  dimethyl  cad¬ 
mium  (DMCd),  Cd(CH3)2>  as  parent  gases  from  which  A1  or  Cd  were  deposited.  These  gases 
have  absorption  bands  with  peaks  near  200  nm  and  measurable  absorptions  for  X  <  260  nm  at 
pressures  of  the  order  of  1  Torr. 


Fig.  II- 9.  Schematic  diagram 
of  the  experimental  system 
used  for  deposition. 


Figure  II- 9  shows  the  experimental  setup  used.  The  output  of  a  CW  argon  ion  laser  operat¬ 
ing  at  514.5  nm  is  frequency  doubled  in  an  ADP  crystal;  the  residual  green  light  is  removed  by 
passing  the  beam  from  the  crystal  through  a  prism  and  the  257.2-nm  radiation  is  focused  on  the 
quartz  exit  window  of  the  3.2-cm-long  stainless  steel  cell  by  either  a  7-cm  focal  length  lens  or 
a  0.5  N.A.  reflecting  microscope  objective.  The  257.2-nm  light  transmitted  by  the  deposit  is 
monitored  by  a  photomultiplier  and  the  resulting  signal  displayed  on  an  X-Y  recorder,  allow¬ 
ing  real-time  monitoring  of  the  growth  of  the  deposit.  Typical  experimental  conditions  were 
~2  x  10~4  W  of  257.2-nm  power,  gas  fills  of  3  Torr  DMCd  or  TMA1  buffered  with  750  Torr  he¬ 
lium,  and  exposure  times  of  1  to  60  sec.  The  focus  of  the  lens  was  located  by  translating  the 
cell  along  the  beam  direction  and  noting  the  position  at  which  the  transmitted  signal  dropped 
most  rapidly.  Deposits  were  examined  by  dark-field  and  Nomarski  microscopy,  by  scanning 
electron  microscope  (SEM),  and  by  Auger  electron  spectroscopy. 

Figure  II- 10  shows  a  Cd  deposit  obtained  by  irradiating  a  3-Torr  DMCd;700-Torr  He  mix¬ 
ture  with  0.10  mW  of  257.2-nm  radiation  for  60  sec  with  the  substrate  out  of  the  focal  plane  of 
the  7-cm  focal  length  lens;  other  deposits  made  with  the  substrate  at  the  focal  plane  produced 
spots  with  lateral  dimensions  as  small  as  3  nm.  This  spot  size  corresponds  to  nearly  diffraction 
limited  performance  by  the  single-element  f/10  lens  and  does  not  indicate  the  potential  resolution 
of  the  process.  Figure  II- 10  shows  that  the  Fresnel  diffraction  pattern  in  the  plane  of  the  win¬ 
dow,  imposed  by  slight  aperturing  of  the  Gaussian  laser  mode,  is  sharply  reproduced  in  the 
metal  deposit.  Examination  of  the  diffraction  rings  obtained  in  this  set  of  depositions  by  SEM 
showed  a  series  of  well -modulated  fringes  with  ring  widths  of  2  (jun.  The  SEM  examination  re¬ 
vealed  the  presence  of  some  secondary  interference  fringes  with  widths  of  0.2  pm;  while  uninten¬ 
tional,  these  fringe  patterns  serve  to  illustrate  that  the  ultimate  resolution  of  the  deposition  pro¬ 
cess  is  in  the  submicrometer  range.  The  SEM  photos  also  showed  that  the  typical  height  of  such 

O 

deposits  was  in  the  1000-  to  2000 -A  range. 
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Fig.  II- 10.  Cadmium  deposit  on  SiC>2  substrate  produced  by  0.1-mW, 

257.2-nm  irradiation  of  3-Torr  DMCd:700-Torr  He  mixture  for  60  sec. 

In  subsequent  experiments  1.0-um  Cd  lines  were  produced  by  translating  the  cell  in  a  direc¬ 
tion  perpendicular  to  the  beam  axis  by  a  motor-driven  micrometer.  In  these  experiments,  the 
reflecting  optics  were  used  and  the  width  of  1.0  ±  0.2  nm  corresponds  closely  to  the  minimum 
spot  size  specified  by  the  manufacturer  of  the  lens,  suggesting  that  smaller  dimensions  will  be 
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possible  with  diffraction  limited  optics.  The  thickness  of  the  deposit  is  estimated  to  be  ~100  A 
from  its  optical  transmission;  greater  thicknesses  are  possible  with  a  slightly  altered  illumina¬ 
tion  geometry.  Writing  speeds  were  limited  by  the  translator  to  ~25  yim/ sec.  A  faster  trans- 
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lator  will  allow  the  UV  laser  intensity  to  be  increased  to  its  maximum  value  (10  greater  than 
at  present)  and  hence  greater  writing  rates  will  be  possible. 

D.  J.  Ehrlich 
T.  F.  Deutsch 
R.  M.  Osgood,  Jr. 


H.  IMAGING  AND  RCS  MEASUREMENTS  OF  SUBMILLIMETER 
MODELED  TACTICAL  TARGETS 


A  submillimeter  system  has  been  developed  for  and  applied  to  imaging  and  radar  cross  sec 
tion  (RCS)  measurements  of  scaled  tactical  targets.  Current  development  of  millimeter  radar 
systems  has  created  a  requirement  for  RCS  measurements  of  targets  over  a  wide  range  of  fre¬ 
quencies.  A  technique  used  for  obtaining  such  data  at  a  specified  frequency  of  interest  vQ  con¬ 
sists  of  measuring  the  power  return  from  a  scale  model  of  the  target  at  a  modeling  frequency 

vm  8iven  by 


where  D  is  a  linear  dimension  of  the  target  and  D  is  the  corresponding  linear  dimension  of 
o  ,,  m 

the  model.  The  modeling  technique  allows  a  large  amount  of  accurate  data  from  complex  tar¬ 
gets,  such  as  a  tank,  to  be  inexpensively  and  rapidly  acquired.  Details  can  be  faithfully  pre¬ 
served  and  model  size  held  to  a  manageable  level  at  scale  factors  between  roughly  20:1  and  100:1 
(Ref.  32). 


Assuming  that  the  range  of  frequencies  over  which  radar  information  is  desired  extends 
from  approximately  10  to  240  GHz,  it  follows  that  the  requisite  frequency  region  for  modeling 
includes  the  entire  submillimeter  (1  mm  to  100  pm)  and  much  of  the  far  infrared  (<100  pm).  A 
particularly  appealing  feature  of  the  modeling  system  described  in  this  report  is  the  use  of  opti¬ 
cally  pumped  lasers  which  provide  hundreds  of  available  frequencies  distributed  throughout  the 
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submillimeter  and  far  infrared.  As  a  result,  a  single  experimental  setup  allows  the  modeling 
of  virtually  any  microwave  frequency  of  interest.  Inexpensive  but  generally  very  accurate  models 
are  often  available  from  commercial  model  manufacturers,  with  suitable  scale  factors  for  these 
measurements. 


MODEL  1 35  1 ) 


Fig.  II- 11.  Schematic  of  the  submillimeter  model  system  showing  the 
monostatic  imaging  configuration. 

The  initial  experiments  described  herein  involved  imaging  studies  of  a  tank.  Figure  II- 11 
shows  the  experimental  setup  of  the  imaging  system.  Submillimeter  radiation  is  directed  by  a 
beam  splitter  through  a  two-dimensional  raster  scanning  system  which  scans  the  beam  over  a 
2 5-cm-dia.f/2  spherical  mirror.  The  near  plane  wave  2-cm-dia.  Gaussian  beam  is  focused  to 
a  waist  and  scanned  over  the  target,  which  is  positioned  in  the  focal  region  of  the  mirror.  Re¬ 
turns  from  the  target  are  detected  by  a  2°K  Ge  bolometer  positioned  behind  the  beam  splitter  so 
as  to  obtain  a  monostatic  view  of  the  target.  The  resulting  output  is  then  displayed  on  a  storage 
oscilloscope. 

Three  methanol  laser  lines  at  1288,  2542,  and  3125  GHz  were  selected  to  model  37,  70,  smd 
90  GHz,  respectively,  thereby  covering  a  substantial  fraction  of  the  millimeter  region  of  interest. 
Figure  II- 12(a)  shows  oscilloscope  traces  of  the  image  of  the  tank  taken  at  a  depression  angle 
of  47°  and  an  aspect  angle  of  0°.  The  outline  of  the  tank  model  as  viewed  by  the  transmitter/ 
receiver  is  shown  in  Fig.  II- 12(b).  Considerable  target  structure  is  evident  in  the  images,  with 
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1  ■ 3125  GHz 
35: 1  SCALE  MODEL 
90  GHz 


(b) 


Fig.  11-12.  (a)  Submillimeter  images  of  the  model  target  at  the  frequencies  and 

orientation  indicated.  The  oscilloscopes  axes  are  driven  by  transducer  voltages 
representing  raster  mirror  positions.  The  returning  signal  is  simply  added  to 
the  vertical  input,  (b)  Outline  of  the  model  target  at  the  indicated  orientation. 
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the  detail  sharpening  with  increasing  frequency,  in  part  as  a  result  of  the  increase  in  resolution. 
The  oscilloscope  traces  demonstrate  that  the  system  signal-to-noise  ratio  is  already  adequate 
for  obtaining  quantitative  measurements  of  the  cross  sections  of  individual  scatterers  on  the 
target.  The  data  shown  are  only  relative  intensities  and  have  not  been  normalized  to  correct 
for  the  variation  in  output  power,  beam-splitter  efficiency,  and  spot  size  with  frequency. 

Determinations  of  absolute  radar  cross  sections  of  individual  scatterers  derived  from  cal¬ 
ibrated  imaging  data  as  well  as  absolute  cross  sections  of  fully  illuminated  targets  are  being 
made  and  compared  with  existing  microwave  data.  The  known  dependence  of  radar  cross  sec¬ 
tions  on  transmitter/receiver  polarization,  resolution  element  size,  and  target  aspect  are  being 
studied  at  frequencies  of  interest  for  millimeter-wave  radar  systems.  Finally,  although  the 
preliminary  measurements  described  herein  were  done  with  CW  sources,  the  existence  of  pulsed, 
optically  pumped  submillimeter  lasers  should  allow  the  extension  of  this  modeling  approach  to 
pulsed  millimeter-wave  radar  systems. 

J.  Wald  man!  P.  E.  Duffy 

H.  R.  Fetterman  T.  G.  Bryant  t 


t  Department  of  Physics,  University  of  Lowell,  Massachusetts, 
t  Group  45. 
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III.  MATERIALS  RESEARCH 


A.  CALCULATED  AND  MEASURED  EFFICIENCIES  OF  SHALLOW- 

HOMOJUNCTION  GaAs  SOLAR  CELLS 

We  have  recently  reported1’^  the  fabrication  of  single -crystal  GaAs  shallow -homojunction 
solar  cells  without  Ga .  A1  As  window  layers  that  have  conversion  efficiencies  of  2i  percent 

1 — X  X  i  i 

at  AMI  (air  mass  1).  These  cells  employ  an  n  /p/p  structure,  prepared  by  chemical  vapor 
deposition  (CVD),  in  which  surface  recombination  losses  are  reduced  because  the  n  -layer  is  so 
thin  that  most  of  the  carriers  are  generated  in  the  p-layer  below  the  junction.  High-efficiency 
cells  have  been  fabricated  on  (100)  substrates  of  both  GaAs  (Ref.  1)  and  Ge  (Ref.  2).  For  the 
same  n+-layer  thickness,  cells  on  both  substrate  materials  have  the  same  conversion  efficiency, 
showing  that  the  solar-cell  properties  of  the  CVD  GaAs  layers  are  the  same  in  both  cases.  The 
substitution  of  Ge  for  GaAs  substrates  will  permit  a  major  decrease  in  Ga  consumption  as  well 
as  a  significant  reduction  in  cell  cost. 

We  have  now  utilized  a  simple  analytical  model  for  the  shallow -homojunction  GaAs  solar 
cells  to  obtain  calculated  results  in  very  good  agreement  with  our  experimental  data  for  the  de¬ 
pendence  of  external  quantum  efficiency  and  conversion  efficiency  on  n-layer  thickness.  The 
efficiency  measurements  were  made  on  cells  with  Ge  substrates,  but  the  results  would  be  the 
same  for  devices  with  GaAs  substrates. 


Fig.III-1.  Schematic  cross  section  of 
GaAs  shallow-homojunction  solar  cell 
on  Ge  substrate. 


i  2 

The  fabrication  of  the  shallow -homojunction  cells  has  been  described  elsewhere.  ’  Figure 
III-l  is  a  schematic  cross  section  of  such  a  cell  on  a  Ge  substrate.  The  p  Ge  substrate  is 
highly  doped  with  Ga  (8  X  iO  cm  ).  The  GaAs  p  buffer  layer  and  p  active  layer  are  doped 

with  Zn  to  5  X  101®  and  1  X  10*7  cm  respectively,  while  the  GaAs  n+-layer  is  doped  with  S 
18  -3 

to  5  X  10  cm  ,  the  highest  S  doping  level  that  has  been  achieved  in  our  CVD  system.  Quan¬ 
tum  and  conversion  efficiency  measurements  were  made  on  four  cells  with  areas  of  about 
2 

0.05  cm  ,  fabricated  side  by  side  on  the  same  wafer.  For  each  cell,  measurements  were  made 

at  several  different  values  of  n"1  -layer  thickness,  which  were  obtained  by  using  a  repeated  anod- 

2 

ization  and  stripping  process  to  accomplish  a  step-wise  reduction  in  thickness.  For  each  mea- 

+  ° 
surement  the  n  -layer  was  anodized  to  form  an  AR  coating  approximately  850  A  thick. 

The  theoretical  model  used  for  our  calculations  assumes  that  one  hole -electron  pair  is  gen¬ 
erated  for  each  solar  photon  absorbed.  The  flow  of  photogene  rated  minority  carriers  across 
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the  p-n  junction  produces  a  photocurrent  consisting  of  three  components.  Some  of  the  holes  gen¬ 
erated  in  the  n+-layer  may  be  lost  by  recombination  at  the  GaAs  surface  as  well  as  by  bulk  re¬ 
combination;  however,  since  the  heterojunction  between  the  heavily  doped  p+  GaAs  and  Ge  does 
not  act  as  a  barrier  to  hole  transport,  all  the  holes  reaching  the  p-n  junction  should  be  collected 
at  the  contact  to  the  Ge  substrate,  thus  contributing  the  first  component  of  the  photocurrent. 
Similarly,  of  the  electrons  generated  in  the  p  and  p+  GaAs  layers,  those  that  reach  the  junction 
are  collected  at  the  contact  to  the  n+-layer  and  contribute  the  second  component  to  the  current. 
Some  of  the  electrons  generated  in  the  p+  GaAs  layer  can  diffuse  into  the  Ge  substrate,  produc¬ 
ing  no  current,  but  electrons  generated  in  the  p  GaAs  layer  are  hindered  from  entering  the  p+ 
region  by  the  back-surface  field  (BSF)  resulting  from  the  difference  in  carrier  concentration 
between  the  p-  and  p+-layers.  The  final  component  of  the  photocurrent  is  due  to  the  minority 
carriers  generated  within  the  junction  depletion  region;  it  is  assumed  that  all  these  carriers  are 
collected  because  of  the  large  drift  field  within  this  region.  Although  some  solar  photons  with 
energies  less  than  the  GaAs  energy  gap  (1.43  eV)  will  be  absorbed  by  the  Ge  substrate,  which 
has  a  much  lower  energy  gap  (0.66  eV),  electrons  generated  in  the  substrate  should  not  contrib¬ 
ute  any  measurable  photocurrent,  since  in  order  to  reach  the  p-n  junction  they  would  have  to 
overcome  a  barrier  of  over  0.7  eV,  the  difference  between  the  energy  gaps  of  GaAs  and  Ge. 

For  the  model  described,  the  external  quantum  efficiency  QEext  as  a  function  of  wavelength 
A  is  given  by 


QEext(*>  = 


J  +  J  +  J  i 
p  n  dr 


(III-l) 


where  Jp  =  photocurrent  density  due  to  holes  generated  in  the  n+-layer,  Jn  =  photocurrent  den¬ 
sity  due  to  electrons  generated  in  the  p-  and  p+-layers,  Jdr  =  photocurrent  density  due  to  carriers 
generated  in  the  depletion  region,  F(A)  =  photon  flux  density  at  A,  and  q  =  electronic  charge.  If 
each  GaAs  layer  is  uniformly  doped,  the  three  photocurrent  densities  are  given  by  the  follow¬ 
ing  equations,  which  are  essentially  the  same  as  those  listed  in  Ref.  3  except  that  different  sets 
of  absorption  coefficients  are  used  for  the  n-  and  p-layers  instead  of  a  single  set  for  both  types. 


qF(l  -  R)  or  L 
J  =  - , — _ G_E 

P  <“nLp-‘' 


>Lp  +  “nV  *  exPl"«nxl  (KpLp  cosh  U~  +  sinh  ir) 


KL  sinh  +  cosh  =2- 

p  p  Lp  s 


-“nLpeXp[_anx> 


qF(l  -  R)  a  L 

J  =  - y— j — e —  exp[-a  (x  +  w)] 

(“p  n  ~  P 


a  L  — 
P  n 


KnLn  (cosh  IT  ~  exPl-«pH'l)  +  sinh  ^  +  apLn  exp[-apH']  j  (III-3) 


KL  sinh  = —  +  cosh  f — 
n  n  L 

n  n 


Jdr  =  qF(1  “  R)  exP(-«nxl  f1  ~  exPt-apw]} 


(III-4) 
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where  R  =  reflectivity  of  the  front  surface;  an,  ap  =  absorption  coefficients  of  n-  and  p-layers, 

respectively;  Lp,  Lp  =  minority-carrier  diffusion  lengths  for  holes  and  electrons,  respectively; 

K  =  S  /D  and  K  =  S  /D  ,  where  S  and  S  are  the  surface  recombination  velocities  for  minor- 
p  p'  p  n  n'  n  p  n 

ity  holes  and  electrons,  respectively,  and  Dp  and  Dn  are  the  minority-carrier  diffusion  coeffi¬ 
cients;  x  =  junction  depth  (essentially  the  thickness  of  the  n+-layer);  w  =  width  of  depletion  re¬ 
gion;  H'  =  H  —  w,  where  H  is  the  effective  thickness  of  the  p-layer.  (Although  F,  R.  <*n,  and 
Op  are  functions  of  X,  in  the  above  equations  and  the  following  discussion  the  functional  symbol 
(X)  has  been  omitted  for  clarity.) 

The  values  of  the  nonadjustable  parameters  appearing  in  Eqs.  (III-2)  to  (III -4)  were  obtained 

4 

as  follows.  The  values  used  for  F  were  those  published  for  AMI.  Values  of  R  were  measured 

O 

on  test  wafers  with  an  anodic  oxide  layer  850  A  thick.  For  <v  and  we  adopted  the  values  re- 

ported3  for  the  absorption  coefficients  of  GaAs  with  n  =  5  x  10  °  cm  and  p  =  1  x  10  cm  , 

respectively.  (For  convenience,  the  same  values  of  op  were  used  for  the  p-  and  p+-layers, 

since  the  actual  coefficients  for  these  two  doping  levels  differ  only  for  X  =  0.8-0. 9  pm,  where 

the  p+  values  are  slightly  lower.)  The  value  used  for  the  electron  diffusion  length  was  L  = 

/  i  _  o  n 

23  pm,  as  reported  for  GaAs  with  p  =  1  X  10  cm”  ;  the  values  calculated  for  QE0xt  are  quite 

insensitive  to  the  exact  value  selected  for  L  ,  as  long  as  the  latter  is  much  greater  than  the 
thickness  of  the  p-layer.  It  was  assumed  that  Kn  =  0,  since  Sn  is  expected  to  be  very  small  be¬ 
cause  of  the  confinement  of  electrons  by  the  BSF.  The  width  of  the  depletion  region  was  calcu- 
7  + 

lated  from  the  carrier  concentrations  in  the  n  -  and  p-layers.  The  depletion  region  lies  almost 
entirely  within  the  p-layer.  Therefore  the  calculated  width  of  the  portion  within  the  p-region, 

0.13  pm,  was  adopted  as  the  value  of  w,  and  the  contribution  of  holes  generated  in  the  n+  deple¬ 
tion  region  to  the  photocurrent  density  is  neglected  in  Eq.  (III-4). 

Since  Lp,  Kp,  and  H  could  not  be  evaluated  independently,  they  were  treated  as  adjustable 
parameters.  A  single  value  for  each  of  these  quantities  was  found  by  performing  computer  cal¬ 
culations  to  achieve  the  best  simultaneous  fit  to  the  six  sets  of  experimental  data  for  ob¬ 

tained  in  measurements  at  six  different  values  of  n+-layer  thickness.  Each  of  these  six  values 

was  found  by  subtracting,  from  the  initial  thickness  of  the  grown  n+-layer,  the  thickness  that 

2 

had  been  removed  by  successive  anodization  and  stripping  steps.  Although  the  thickness  re¬ 
moved  was  determined  quite  accurately  from  the  anodization  voltages  used,  the  initial  n+-layer 
thickness  was  less  certain.  Therefore  this  initial  thickness  was  also  treated  as  an  adjustable 

parameter.  The  initial  value  giving  the  best  fit  to  all  six  sets  of  QE  .  data  was  2150  A,  100  A 

2  exl  + 
greater  than  the  value  estimated  from  the  anodization  results.  The  values  of  n  -layer  thick¬ 
ness  used  throughout  the  remainder  of  this  section  were  obtained  by  using  this  best-fit  initial 
thickness. 

In  Fig.III-2  the  calculated  curves  for  QEgxt  as  a  function  of  X  are  compared  with  the  ex¬ 
perimental  points  for  four  values  of  x.  The  overall  agreement  is  seen  to  be  excellent,  although 

O 

for  x  =  450  A  at  the  shorter  wavelengths  the  calculated  quantum  efficiency  values  are  somewhat 

O 

higher  than  the  measured  ones.  Good  agreement  was  obtained  for  x  =  300  and  1300  A,  the  re- 

suits  for  which  have  been  omitted  from  Fig.  III-2  for  clarity,  except  that  the  calculated  QE„  . 

o  exx 

values  for  x  =  300  A  are  significantly  higher  than  the  measured  ones  for  X  <  0.6  pm. 

The  following  adjustable  parameter  values  were  used  in  calculating  the  QE  ,  curves  in 

ext  7  -1 

order  to  give  the  best  fit  to  experiment:  H  =  4  pm,  Lp  =  0.05  pm,  and  Kp  =  1.7  X  10  cm  . 

An  even  higher  value  for  H,  the  effective  p-layer  thickness,  would  give  a  better  fit  to  the  QEejrt 
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Fig.  III-2.  External  quantum  efficiency 
as  a  function  of  wavelength  for  GaAs 
shallow -homojunction  solar  cells  with 
n+-layer  thickness  of  450,  750,  1000. 

O 

and  1550  A.  The  solid  curves  were  cal¬ 
culated,  while  the  points  were  measured. 


data  for  0.8  to  0.85  pm,  but  such  a  value  would  be  meaningless  since  4  pm  is  the  total  thickness 
of  the  p-  and  p+-layers.  A  value  closer  to  the  p-layer  thickness  (2  pm)  would  be  expected  be¬ 
cause  not  all  the  electrons  generated  in  the  p+-layer  should  reach  the  p-n  junction.  This  dis¬ 
crepancy  suggests  that  the  actual  absorption  coefficients  for  the  p-layer  may  be  larger  than  the 
literature  values,  since  the  calculations  would  yield  a  smaller  value  of  H  if  larger  values  were 
used  for  a  . 

P  j.  18-3 

The  Lp  value  for  the  n  -layer,  which  was  doped  with  S  to  5  X  10  cm  ,  is  considerably 

lower  than  the  values  of  0.1  and  0.3  pm  reported  for  GaAs  layers  doped  to  the  same  level  with 

Te  (Ref.  8)  or  Sn  (Ref.  9),  respectively.  The  reduction  in  Lp  could  be  explained  if  precipitates 

or  defects  were  present  in  the  n+-layer  as  a  result  of  doping  the  layer  with  S  to  a  level  near  the 

solubility  limit.  This  effect  could  also  explain  the  fact  that  the  electron  Hall  mobility  mea- 
+  2-1-1 

sured  for  the  n  -layer  is  only  about  1000  cm  V  sec  ,  compared  to  the  value  of  about 
2-1-1  8 

2000  cm  V  sec  reported  for  layers  with  the  same  level  of  Te  doping.  If  it  is  assumed  that 

the  drift  mobility  of  electrons  is  equal  to  the  Hall  mobility,  that  the  ratio  of  hole  mobility  to 

electron  mobility  in  the  n+-layer  is  1/20  [a  reasonable  value  for  GaAs  (Ref.  7)1.  and  that  the 

Einstein  relation  is  applicable,  then  Dp  =  (kT/q)pp  =  1.3  cm  sec  .  Combining  this  Dp  value 

with  the  value  obtained  for  K  =  S  /D  gives  an  estimated  surface  recombination  velocity 
7  P  P  P  ,  7-1 

S  »  2  X  10  cm  sec  for  the  anodized  n+-layer,  consistent  with  the  value  of  10  cm  sec”  re- 
P  in 

ported  for  a  bare  n  GaAs  surface. 


By  using  the  same  parameter  values  used  in  the  QEext  calculations,  we  have  calculated  the 
maximum  conversion  efficiency  V  max  at  AMI  as  a  function  of  n+-layer  thickness  for  cells  with 
the  structure  shown  in  Fig.  Ill— 1.  The  equations  employed  are  given  in  Refs.  3  and  11.  For  the 
values  of  dark  saturation  current  density  Jq  appearing  in  these  equations,  we  used  calculated 


values  of  the  injected  current  density  J.^.  For  our  cells,  J.^  should  closely  approximate  JQ, 
since  the  measured  diode  factors  are  close  to  unity.  The  values  calculated  for  J.  .  depend 
strongly  on  the  value  adopted  for  the  intrinsic  carrier  concentration  m.  In  our  calculations  we 
used  nt  =  1.8  X  106  cm”3  at  300  K  (Refs.  12  and  13). 

The  dependence  of  AMI  conversion  efficiency  on  n+-layer  thickness  is  shown  in  Fig.  III-3, 


where  the  measured  values  of  9exp  (represented  by  points)  for  four  different  n+/p/p+  cells  are 
compared  with  the  calculated  values  of  9  max  (represented  by  the  solid  line).  With  decreasing 
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Fig.  III-3.  Conversion  efficiency  at  AMI 
as  a  function  of  n+-layer  thickness  for 
GaAs  shallow -homojunction  solar  cells. 
The  measured  values  are  represented  by 
points,  the  calculated  values  of  maxi¬ 
mum  conversion  efficiency  r/  max  and  of 

0.9  n  by  solid  and  dashed  lines, 

'max  J 

respectively. 


thickness,  n  increases  from  17  percent  at  1600  A  to  24  percent  at  300  A.  Even  at  the  latter 

thickness,  about  7  percent  of  the  total  photocurrent  is  due  to  the  collection  of  holes  generated  in 

the  n+-layer.  The  n  calculations  do  not  take  account  of  losses  resulting  from  such  effects 
max 

as  contact-finger  shadowing  and  reduction  in  fill  factor  due  to  series  resistance.  Because  of 
such  losses,  we  expect  the  tj  values  for  our  cells  to  approximate  0.9  r>,„  (represented  by 
the  dashed  line  in  Fig.III-3).  Within  the  limits  of  experimental  error  this  is  generally  the  case, 
although  the  V  cxp  values  do  not  appear  to  decrease  as  rapidly  as  V  max  with  increasing  n+-layer 
thickness. 


In  summary,  by  means  of  calculations  based  on  a  simple  model  we  have  obtained  good  fits 
to  the  measured  values  for  the  external  quantum  efficiency  and  conversion  efficiency  of  GaAs 
solar  cells  with  the  shallow-homojunction  n+/p/p+  structure.  These  calculations  have  yielded 
values  for  several  material  properties  of  the  GaAs  layers  composing  these  cells.  Having  dem¬ 
onstrated  the  applicability  of  the  model,  we  can  now  use  it  in  optimizing  cell  design.  Preliminary 
calculations  indicate  that  still  higher  conversion  efficiencies  would  be  obtained  for  n+-layers  with 
higher  minority-carrier  diffusion  lengths  and  lower  surface  recombination  velocities. 


J.  C.  C.  Fan 
C. O.  Bozler 
B.  J.  Palm 


B.  CALCULATION  OF  CRYSTALLIZATION  FRONT  VELOCITY  FOR  AMORPHOUS 

SEMICONDUCTOR  FILMS  SCANNED  WITH  A  CW  LASER 

In  experiments  on  crystallization  of  amorphous  semiconductor  films  by  scanning  with  the 
slit  image  of  a  CW  laser,  we  have  observed  a  number  of  unusual  phenomena,  including  the  for¬ 
mation  of  periodic  structural  features  on  the  film  surface,  pulsations  of  the  film  temperat  tr.' 
during  scanning,  and  runaway  crystallization  of  the  entire  film  following  first  contact  with  the 
laser  image.  By  means  of  preliminary  computer  calculations  of  the  crystallization  front  veloc¬ 
ity,  we  have  shown  that  all  these  phenomena  can  be  attributed  to  the  liberation  of  heat  accompany¬ 
ing  the  transformation  from  the  amorphous  to  the  crystalline  state. 
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Although  semiconductor  films  deposited  in  the  amorphous  state  will  remain  in  this  state 
indefinitely  at  room  temperature,  they  are  metastable  with  respect  to  the  crystalline  form.  When 
the  film  temperature  is  raised  sufficiently,  transformation  to  the  crystalline  form  takes  place, 
accompanied  by  the  liberation  of  heat.  The  transformation  rate  increases  exponentially  with  ab¬ 
solute  temperature,  so  that  over  a  narrow  temperature  interval  the  rate  changes  by  orders  of 
magnitude.  For  a  zeroth-order  calculation  of  the  motion  of  the  amorphous -crystalline  phase 
boundary  produced  by  laser  scanning,  we  can  therefore  assume  that  the  transformation  occurs 
instantaneously  when  the  amorphous  film  reaches  a  sharply  defined  critical  temperature  T  .  We 
also  neglect  the  temperature  dependence  of  the  optical  and  thermal  properties  of  the  film,  as  well 
as  the  difference  in  these  properties  between  the  amorphous  and  crystalline  states.  With  these 
simplifying  assumptions,  an  integral  equation  formalism14  can  be  used  to  describe  the  motion  of 
the  phase  boundary. 

The  geometry  assumed  for  the  laser  crystallization  calculation  is  shown  schematically  in 
Fig.  Ill— 4.  The  semiconductor  film,  which  is  deposited  on  a  thick  substrate,  is  of  infinite  ex¬ 
tent  in  the  y  and  z  directions  and  so  thin  that  its  temperature  is  constant  in  the  x  direction. 


Fig.  111-4.  Schematic  diagram  of  geometry 
assumed  for  theoretical  description  of  a 
laser  crystallization  experiment. 


The  laser  slit  image  is  of  infinite  length  in  the  z  direction  and  moves  at  a  velocity  v  from  left 
to  right  in  the  y  direction.  At  t  =  0,  the  phase  boundary  is  located  at  y  ,  with  the  crystalline 
phase  to  the  left  (y  <  yQ)  and  the  untransformed  amorphous  phase  to  the  right  (y  >  yQ).  Prior  to 
t  =  0,  the  laser  image  is  moving  through  the  region  y  <  y  ,  carrying  with  it  a  steady-state  tem¬ 
perature  profile  (y  —  vt).  At  t  =  0,  the  temperature  at  the  phase  boundary  reaches  Tc  and  the 
boundary  begins  to  move  toward  the  right,  with  heat  being  liberated  at  a  rate  per  unit  volume  of 
LpY(t),  where  L  is  the  latent  heat  of  crystallization  of  the  semiconductor,  p  is  its  density,  and 
Y(t)  is  the  position  of  the  boundary.  As  a  further  simplification,  the  heat  diffusing  into  the  sub¬ 
strate  and  the  heat  lost  by  radiation  from  the  film  surface  are  neglected.  Then  the  temperature 
T(y,  t)  at  any  point  along  the  film  at  time  t  is  given  by  the  integral  relation 


T(y,t)  =  Tj(y  —  vt)  +  i 


Y(t') 
(t  - 


d" 


(III-5) 
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where  C  is  the  specific  heat  of  the  film,  it  =  K/Cp  defines  its  thermal  diffusivity,  and  K  is  its 
thermal  conductivity.  An  integral  equation  for  Y(t)  can  be  obtained  by  using  the  condition  that 
the  temperature  at  the  phase  boundary  is  T  ,  or 

T[Y(t),  t]  =  Tc  .  (III-6) 

For  purposes  of  calculation,  it  is  convenient  to  rewrite  Eq.  (III-5)  in  the  frame  of  reference 
moving  with  the  laser  image.  For  this  purpose  we  introduce  the  position  variable  u(t)  =  y(t)  —  vt, 
where  u(t)  is  measured  from  the  center  of  the  laser  image  as  origin.  The  temperature  Tj  is 
modeled  in  the  form 


T1(u)  =  Tfa  +  Tl  e 


-(u/a) 


(III -7 ) 


where  is  a  uniform,  time -independent  background  temperature,  and  the  temperature  contri¬ 
bution  due  to  the  laser  is  described  by  a  Gaussian  of  width  a. 

Finally,  introducing  normalized  quantities  we  can  write  the  integral  equation  for  the  motion 
of  the  phase  boundary  as 


i  =  a  exp{-[S(T)]2}  +  77  C  - 4  J2  exp 
Jo  (t-t  'V' 


—  [S (r )  -S(T')  +  V (t  -  T' 


where  a  =  Tj  /(Tc  -  Tb),  7)  =  L/Ctt1/2(Tc  -  Tb),  r  =  4*t/a2,  V  £  av/4ic,  S  =  U/a,  S(r')  = 
dS(T')/dr'  =  (l/a)  (dU(t')/dt']  (AV /At  ’),  and  U  is  the  position  of  the  phase  boundary  measured 
from  the  center  of  the  laser  image. 

Equation  (III -8 )  has  been  solved  numerically  by  a  preliminary  method  to  obtain  S  as  a 
function  of  t  for  representative  values  of  a,  77,  and  V,  with  the  boundary  conditions  that 
S(t)  =  0  at  t  =  0  and  that  S(0)  is  given  by  1  =  a  exp  {  — [S(0)]2} .  Figure  I II - 5  shows  plots  of  S  vs 
r  for  o’  =  1.5,  V  =  1.52,  and  three  values  of  77.  For  each  value  of  77 ,  S  initially  increases  rap¬ 
idly  (i.e.,  the  amorphous -crystalline  phase  boundary  moves  quickly  ahead  of  the  advancing  laser 
image)  because  the  heat  liberated  by  the  phase  transformation  raises  the  temperature  ahead  of 
the  boundary,  accelerating  the  forward  motion  of  the  boundary.  As  the  boundary  moves  away 
from  the  laser  image,  the  contribution  of  the  laser  to  the  temperature  ahead  of  the  boundary 
decreases  rapidly.  For  77  =  0.19,  the  heat  liberated  by  the  transformation  is  then  insufficient 
to  keep  the  boundary  velocity  greater  than  the  laser  velocity,  and  S  therefore  decreases.  With 
the  resulting  approach  of  the  laser  image  toward  the  boundary,  the  laser  contribution  to  the 
boundary  temperature  increases,  eventually  becoming  sufficient  to  once  more  accelerate  the 
motion  of  the  boundary  and  causing  S  to  increase  again.  This  cycle  is  repeated  indefinitely, 
resulting  in  the  oscillations  in  S  seen  for  77  =  0.19  in  Fig.III-5. 

As  7)  increases,  the  heat  liberated  by  the  phase  transformation  increases  relative  to  the 


heat  required  to  raise  the  film  temperature  from  the  background  temperature  Tb  to  the  critical 
temperature  Tc<  Therefore,  the  initial  acceleration  of  the  phase  boundary  is  greater,  and  the 
boundary  moves  farther  beyond  the  laser  image  before  it  decelerates.  This  trend  is  illustrated 
by  the  initial  portion  of  the  calculated  curve  for  77  =  0.35  in  Fig.  III-5.  The  remainder  of  this 
curve  has  been  dotted  because  the  solution  obtained  beyond  t  =  0.085  is  nonphysical.  Thus,  the 
decrease  in  the  calculated  value  of  S  for  0.085  <  t  <  0.15  corresponds  to  a  backward  motion  of 
the  phase  boundary,  which  would  imply  reconversion  of  the  crystalline  phase  to  the  amorphous 
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Fig.III-5.  Calculated  curves  of  S(r)  vs  t 
for  parameter  values  a  =  1.5,  V  =  1.52, 
with  n  =  0.19,  0.35.  and  0.8. 


state  and  a  consequent  reabsorption  of  the  latent  heat.  In  reality  the  boundary  comes  to  rest  at 
t  ~  0.085  and  then  remains  fixed  until  the  laser  image  approaches  closely  enough  to  again  initiate 
forward  motion.  This  behavior  is  indicated  schematically  by  the  dashed  curve  in  Fig.  111-5,  the 
linear  portion  of  which  covers  the  time  interval  when  the  boundary  is  fixed  and  therefore  has  the 
slope  |  S |  =  V.  (It  should  be  noted  that  during  this  time  interval  the  boundary  temperature  ini¬ 
tially  drops  below  Tc  and  then  increases  to  T  ,  at  which  time  the  boundary  begins  to  move  for¬ 
ward.  For  the  motion  represented  by  the  solid  lines  in  Fig.  III-5,  the  temperature  at  the 
boundary  is  always  Tc.)  The  physically  correct  solution  for  the  boundary  motion  can  be  ob¬ 
tained  by  solving  Eq.  (III-8)  with  the  additional  constraint  that  whenever  [S(t')  +  V)  <  0,  this 
quantity  should  be  set  equal  to  zero. 

When  7)  becomes  large  enough,  the  heat  liberated  during  crystallization  is  sufficient  to 
sustain  the  transformation,  causing  the  boundary  to  "run  away"  from  the  laser  image.  This 
situation,  which  is  illustrated  by  the  curve  for  ij  =  0.8  in  Fig.III-5,  corresponds  to  the  obser¬ 
vation  that  the  entire  film  is  crystallized  following  first  contact  with  the  laser  image. 

The  zeroth-order  model  described  does  not  predict  observable  structural  changes  in  the 
laser-treated  films,  since  optical  and  thermal  properties  are  assumed  identical  for  the  amor¬ 
phous  and  crystalline  phases.  Nevertheless,  the  qualitative  features  of  laser  crystallization 
experiments  can  be  understood  on  the  basis  of  the  model  calculation  if  the  more  realistic  prop¬ 
erties  of  the  films  are  taken  into  account  in  the  next  approximation.  Among  the  most  important 
considerations  is  the  fact  that  the  rate  of  conversion  of  amorphous  to  crystalline  material  is  a 
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smooth  function  of  temperature,  rather  than  occurring  abruptly  at  T  .  The  actual  morphology 
of  different  regions  of  a  laser  -treated  film  can  be  expected  to  depend  on  their  rates  of  transfor¬ 
mation  and  therefore  on  their  thermal  history.  Furthermore,  crystalline  material  absorbs  less 
laser  energy  than  amorphous  material,  so  that  even  though  all  the  material  is  exposed  directly 
to  the  laser  radiation  after  the  phase  boundary  has  passed  by,  this  direct  exposure  may  not  fur¬ 
ther  modify  the  morphology  of  the  film.  This  suggests  that  the  periodic  structural  features  ob¬ 
served  on  laser-crystallized  films  can  be  attributed  to  periodic  fluctuations  in  the  velocity  of 
the  phase  boundary  like  those  shown  by  the  curves  for  ij  =  0.19  and  0.35  in  Fig.  Ill— 5.  If  so, 
comparison  of  these  curves  indicates  that  the  features  should  be  more  widely  spaced  and  more 
prominent  for  higher  values  of  rj  (i.e.,  higher  T^).  Preliminary  observations  on  laser- 
crystallized  films  of  Ge  on  fused  silica  appear  to  be  consistent  with  this  prediction.  Further¬ 
more,  for  sufficiently  high  values  of  t),  the  calculated  phase  boundary  alternately  moves  rapidly 
and  comes  to  rest,  resulting  in  large  fluctuations  in  the  rate  of  heat  liberation  and  therefore  in 
temperature.  This  could  explain  the  observation  of  periodic  fluctuations  in  the  color  tempera¬ 
ture  of  films  during  laser  crystallization. 


Fig.  I II -6 .  Microphotograph  of  periodic  structure  observed  after  laser 
crystallization  of  a  film  of  Ge  on  a  fused  silica  substrate. 

Figure  III-6  is  a  photomicrograph  showing  structural  features  with  a  spatial  period  of 
~300  p.m  produced  during  laser  crystallization  of  an  amorphous  Ge  film  on  a  fused  silica  sub¬ 
strate.  Assuming  that  such  features  are  due  to  oscillations  in  phase  boundary  velocity  like 
those  shown  in  the  curves  of  Fig.lII-6,  their  spatial  period  AY  can  be  estimated  from  these 
curves  by  using  the  relationship  iY/a  =  VAt,  where  at  is  the  temporal  period  of  the  velocity 
oscillation.  For  n  =  0.19  and  0.35,  the  values  of  AY/a  are  0.1  and  0.4,  respectively.  A  rough 
estimate  of  the  parameter  a  for  Ge  on  fused  silica  is  10  1  cm,  yielding  AY  values  of  100  and 
400  pm,  respectively.  In  view  of  the  preliminary  nature  of  the  calculations,  this  degree  of 
agreement  with  the  experimental  observations  is  encouraging. 

H.J.  Zeiger  J.C.C.  Fan 

B.  J.  Palm  R.  L.  Chapman 
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IV.  MICROELECTRONICS 


A.  CRYSTALLOGRAPHIC  ORIENTATION  OF  SILICON  ON  AMORPHOUS  SUBSTRATE 
USING  ARTIFICIAL  SURFACE-RELIEF  GRATING  AND  LASER  CRYSTALLIZATION 

In  earlier  work,  surface -relief  gratings  with  square-wave  cross  section  in  amorphous  SiO- 

1  * 

were  used  to  induce  oriented  growth  of  KC1  crystallites  grown  from  a  water  solution,  and 

2-4 

to  align  overlayer  films  of  nematic  and  smectic  A  liquid  crystals.  In  that  work  and  other 
papers5,  ^  the  authors  speculated  that  control  of  overlayer  crystallographic  orientation  by  means 
of  artificial  surface-relief  structure  should  be  applicable  to  a  broad  range  of  overlayer/substrate 
combinations  and  deposition  methods.  Here  we  report  uniform  crystallographic  orientation  of 
thin  silicon  films  on  amorphous  fused  silica  substrates,  achieved  by  laser  crystallization  of 
amorphous  silicon  deposited  over  a  surface-relief  grating  in  the  substrate.  We  propose  that 

such  examples  of  orientation  induced  by  artificial  surface  patterns  be  called  "graphoepitaxy."* 

3 

A  theoretical  model  predicts  that  film  formation  methods  that  yield  textured  polycrystalline 
films  (i.e.,  films  in  which  the  individual  crystal  grains  have  one  particular  plane  parallel  to  the 
substrate  surface  and  random  orientation  otherwise)  on  smooth  amorphous  substrates  should 
yield  uniformly  oriented  films  if  such  film  formation  is  carried  out  over  an  appropriate  surface - 
relief  structure  in  the  amorphous  substrate.  For  example,  polycrystalline  films  with  (100)  tex¬ 
ture  [i.e.,  with  grains  having  (100)  planes  parallel  to  the  substrate  surface  but  random  orienta¬ 
tion  in  the  surface  plane)  should  become  uniformly  oriented  on  a  surface-relief  grating  with 
square -wave  cross  section  if  the  spatial  period  of  the  grating  is  small  compared  with  the  normal 
grain  size.  Specifically,  the  model,  which  is  based  on  equilibrium  thermodynamics,  states  that 
such  uniform  orientation  is  a  configuration  of  minimum  free  energy.  Polycrystalline  films  of 
(111)  texture,  on  the  other  hand,  would  require  a  relief  structure  having  facets  that  intersect 

at  70.5°  and/or  109.5°.  Relief  structures  with  square-wave  cross  section  are  more  easily  fab- 

3  7  8 

ricated  using  current  techniques,  '  ’  and  thus  one  is  led  to  seek  film  formation  methods  that 

yield  (100)  or  (110)  textured  silicon.  It  is  known  that  laser  crystallization  of  silicon  films  yields 

9  10 

grains  several  micrometers  in  diameter.  ’  We  found  that  under  certain  conditions  laser  crys¬ 
tallization  of  amorphous  silicon  over  fused  silica  yields  a  nearly  perfect  (100)  texture. 

Several  experiments  were  performed  to  determine  the  texture  of  laser-crystallized 
amorphous  and  polycrystalline  silicon,  deposited  by  a  variety  of  methods,  on  optically  polished 
fused  silica  substrates.  In  some  cases  the  polished  fused  silica  was  coated  with  SiO^  deposited 
by  a  CVD  process.  With  amorphous  silicon  films  500-nm  thick,  deposited  at  610°C  in  a  com- 
t  N2 

mercial  CVD  reactor*  using  the  reaction  SiH^  — »  Si  +  we  were  able  to  obtain  100  percent 

of  a  (100)  texture  after  laser  crystallization,  as  determined  by  x-ray  diffractometry.  The 
maximum  deviation  of  (100)  from  the  substrate  plane  was  about  3°,  as  determined  by  reflection 
electron  diffraction  (RED).  With  silicon  deposited  by  thermal  evaporation  and  by  the  reaction 
SiH^  -*  Si  +  ZH^  in  a  low-pressure  CVD  reactor*  we  observed  100  percent  (111)  and  relatively 
weak  (100)  textures,  respectively.  The  laser  crystallization  was  done  in  air  at  room  tempera¬ 
ture  using  an  argon-ion  laser,  operating  multiline  at  a  power  level  between  6  and  7  W,  focused 
with  a  60-mm  focal  length  lens  placed  45  mm  from  the  silicon  film.  For  a  single  linear  scan. 


t  From  the  Greek:  "grapho"  -  to  write  or  incise;  "epi"  -  upon;  "taxis"  -  arrangement,  order. 
{General  Signal/Tempress,  Sunnyvale,  California. 
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this  beam  produced  a  band  of  crystallization  about  380  pm  wide.  Samples  were  scanned  past 
the  beam  at  10  mm/sec  or  slower  in  a  raster  fashion,  with  successive  horizontal  scans  stepped 
vertically  by  12  pm.  After  multiple  raster  scans  of  a  silicon  film  in  the  laser  beam,  poly¬ 
crystalline  grains  grew  to  10  to  100  pm  diameter,  indicating  that  uniform  crystollographic 
orientation  should  be  obtainable  with  square-wave  gratings  having  spatial  periods  of  several 
micrometers. 

To  obtain  square -wave  gratings  in  fused  silica  substrates  we  first  coated  the  substrates 
with  a  film  of  chromium  20  nm  thick.  Using  conformable  photomask  lithography,11  a  3.8-pm 
spatial-period  grating  was  then  exposed  and  developed  in  AZ  1350B  photoresist,*  and  the  chro¬ 
mium  etched  in  a  chemical  etchant.  After  removal  of  the  photoresist,  a  square-wave  grating 
was  etched  into  the  fused  silica  by  reactive-ion  etching  in  CHF^  gas  using  the  chromium  grating 
as  a  mask.  '  '  Finally,  the  chromium  was  removed  in  a  chemical  etchant.  This  procedure 
yields  square-wave  gratings  with  flat  tops,  flat  bottoms,  and  sidewalls  that  deviate  a  maximum 
of  about  6°  from  the  vertical.  The  radii  of  curvature  at  the  cornex's  where  the  sidewalls  join 
the  *ops  and  groove  bottoms  are  about  5  nm.  The  deviation  of  a  sidewall  from  an  ideal  planar 
facet  depends  primarily  on  the  edge  ripple  in  the  original  photomask  and  any  additional  smooth¬ 
ing  or  roughening  that  occurs  during  the  chemical  etching  of  the  chromium.  Typically  the  de¬ 
viation  is  about  ±100  nm.  In  the  experiments  reported  here  groove  depths  were  19,  60,  and 
100  nm.  Substrates  were  cleaned  just  prior  to  depositing  silicon  on  them,  the  final  cleaning 
step  being  irradiation  for  about  40  min.  in  air  with  a  far-UV  lamp.12  The  lamp  generates  ozone 
and  is  very  effective  in  removing  any  residual  organic  contamination. 

For  the  crystallographic  orientation  experiments,  about  500  nm  of  amorphous  silicon  was 
deposited  by  the  CVD  process,*  referred  to  above,  over  fused  silica  substrates  1.5  mm  thick 
and  38  mm  in  diameter.  The  3.8-pm  spatial-period  grating  was  100  nm  deep  and  covered  a 
square  area  about  15  mm  on  a  side.  The  silicon  was  laser  crystallized,  both  inside  and  outside 
the  grating  area  (using  the  scanning  procedure  described  above),  at  laser  powers  of  6  and  7  W. 
After  a  single  raster  scan  at  6  W,  the  amorphous  silicon  became  polycrystalline,  with  (100) 
texture  and  grains  10  to  20  pm  in  diameter,  but  no  orientation  relative  to  the  grating  was  mea¬ 
surable.  After  a  second  raster  scan,  grains  had  grown  to  about  100  pm  dia.,  and  there  was  a 
definite  preferred  orientation  relative  to  the  grating  as  determined  by  RED.  After  four  raster 
scans,  a  nearly  perfect  orientative  relation  to  the  grating  was  obtained,  where  the  <100>  direc¬ 
tions  in  the  silicon  are  parallel  to  the  grating  and  perpendicular  to  the  substrate  surface,  as 
illustrated  schematically  in  Fig.  IV-l(b).  To  verify  that  the  grating  geometry  (and  not  the  direc¬ 
tion  of  the  laser  scan)  determines  the  silicon  orientation,  scans  were  made  parallel,  perpendic¬ 
ular,  and  at  45°  to  the  grating.  Orientation  was  found  to  be  independent  of  laser  scan  direction. 
As  depicted  in  Fig.  IV-l(a),  the  deposit  of  amorphous  silicon  conforms  to  the  contour  of  the 
underlying  grating.  With  laser  crystallization,  however,  there  is  substantial  redistribution  of 
material;  the  top  surface  of  the  silicon  becomes  roughened  but  it  is  generally  planar  with  no 
evidence  of  even  a  shallow  grating.  Outside  the  grating  area  the  film  consists  of  single 
crystal  grains,  10  to  100  pm  in  diameter,  that  have  a  (100)  texture,  with  random  orientation 
in  the  plane  of  the  substrate.  Figure  IV-2  shows  a  comparison  of  the  RED  patterns  obtained 
inside  and  outside  the  grating  area.  The  small  arcs  in  Fig.  IV-2  indicate  a  slight  deviation 


t  Shipley  Company,  Newton,  Massachusetts. 

1  General  Signal/Tempress,  Sunnyvale,  California. 
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Fig.  IV-1.  Schematic  illustration  of  (a)  amor¬ 
phous  silicon  film  deposited  by  a  CVD  process 
over  a  surface -relief  grating  in  fused  silica  and 
(b)  uniformly  oriented  film  of  silicon  obtained 
by  laser  crystallization  of  amorphous  silicon 
over  relief  structure.  Note  that  in  (b)  the  top 
surface  of  the  silicon  does  not  follow  the  con¬ 
tour  of  the  grating,  as  in  (a).  In  practice,  the 
surface  of  the  silicon  in  (b)  has  a  roughness  of 
several  10  nm. 


of  (100)  texture  from  one  area  in  the  film  to  another.  Such  deviation  is  present  both  inside  and 
outside  the  grating  area. 

Figure  IV-3  shows  a  transmission  diffraction  pattern  taken  in  a  transmission  electron 
microscope  of  a  silicon  film  laser  crystallized  over  the  grating.  The  silicon  was  separated 
from  the  fused  silica  substrate  by  immersion  in  a  HF  solution.  The  Kikuchi  lines  in  the  dif¬ 
fraction  pattern  indicate  a  highly  ordered  crystalline  structure.13 

If  the  laser  power  is  increased  to  7  W,  grains  having  a  variety  of  shapes  up  to  200  pm  dia. 
form  on  the  first  raster  scan,  but  no  crystallographic  orientation  is  observed  in  the  grating 
area.  This  absence  of  orientation  persists  despite  repeated  scans.  Removal  of  the  silicon  in 
hydrazine  indicates  some  damage  to  the  grating  at  the  7-W  power  level.  Also,  at  7  W  the 
silicon  appears  to  be  melted  by  the  laser  irradiation. 

The  cross-sectional  profile  of  the  grating  plays  an  important  role  in  controlling  crystallo¬ 
graphic  orientation.  If  the  depth  of  the  grating  is  reduced  to  60  nm,  the  crystallographic  orien¬ 
tation  degrades  somewhat,  and  at  a  depth  of  19  nm  no  evidence  of  orientation  relative  to  the 
grating  is  observed  by  RED.  If  the  bottom  corners  of  the  grating  are  rounded  to  a  radius  of 
curvature  of  about  80  nm,  the  orientation  effect  of  the  grating  is  substantially  degraded.  This 
is  understandable  qualitatively  from  the  theoretical  model,  and  is  in  agreement  with  similar 
results  obtained  with  KC1  where  orientation  would  not  occur  on  gratings  with  rounded  corners.3 

We  have  demonstrated  uniform  crystallographic  orientation  of  silicon  films  on  amorphous 
fused  silica  substrates  achieved  by  laser  crystallization  over  a  surface-relief  grating  with 
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Fig.  IV-2.  (a)  RED  pattern  showing  the  (100)  texture  of  laser-crystallized 

silicon  film  on  smooth  amorphous  fused  silica  substrate.  Lengths  of  the 
arcs  in  the  pattern  indicate  the  maximum  deviation  of  the  (100)  planes  from 
parallel  with  the  substrate  is  about  3°.  Pattern  is  independent  of  rotation 
about  the  substrate  normal,  indicating  absence  of  any  in-plane  orientation, 
(b)  RED  pattern  of  silicon  film  laser  crystallized  over  a  3.8-pm  spatial  pe¬ 
riod,  100  nm  deep,  square-wave  grating  in  amorphous  fused  silica.  Elec¬ 
tron  beam  is  parallel  to  the  grating  lines.  Observed  (100)  pattern  indicates 
uniformly  oriented  silicon  having  <100>  directions  parallel  to  the  grating 
lines  and  perpendicular  to  the  substrate  surface,  (c)  RED  pattern  of  same 
sample  as  in  (b)  but  with  electron  beam  at  45°  to  the  grating  lines.  Diffrac¬ 
tion  pattern  is  (110).  Appearance  of  arcs  rather  than  spots  in  (b)  and  (c) 
indicates  the  film  consists  of  single-crystal  regions  where  the  (100)  plane 
deviates  a  maximum  of  about  3°  with  respect  to  the  substrate  surface. 
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Fig.  IV-3.  Transmission  electron  diffraction  pattern  of  500-nm-thick 
silicon  film  laser  crystallized  over  a  3.8-pm  period  square -wave  grat¬ 
ing  in  amorphous  fused  silica.  Electron  beam  was  perpendicular  to 
silicon  film  yielding  a  (100)  diffraction  pattern.  Kikuchi  lines  visible 
in  the  pattern  indicate  a  high  degree  of  crystalline  perfection. 


square-wave  cross  section.  This  opens  the  possibility  of  new  methods  of  preparing  films  for 

microelectronic  devices  and  solar  cells.  ^  ^  Qeis 

D.  C.  Flanders 
H.  I.  Smith 


B.  CHARGE-COUPLED  DEVICES:  SAW/CCD  PROGRAMMABLE  MATCHED  FILTER 

We  report  here  a  new  device,  the  SAW/CCD  programmable  matched  filter,  which  is  capable 
of  correlating  an  analog  SAW  signal  of  up  to  40-MHz  bandwidth  and  3.5-psec  duration  against 
300  samples  of  an  analog  reference  previously  loaded  into  a  CCD  shift  register.  The  program¬ 
mable  matched  filter  is  the  third  analog  signal  processing  function  to  be  demonstrated  by  this 

new  class  of  integrated  SAW/CCD  devices.  By  adding  the  programmable  matched  filter  to  the 

14  15 

previously  reported  buffer  memory  and  accumulating  correlator,  we  have  a  family  of  unique 
devices  which  combines  the  wide  bandwidth  of  the  SAW  with  the  long  and  flexible  time-handling 
capability  of  the  CCD,  and  thereby  promises  significant  advantages  over  other  technologies  in 
many  signal  processing  applications. 

The  construction  details  for  this  device  are  similar  to  those  reported  earlier  for  the  SAW/ 
CCD  fast-in,  slow-out  buffer  memory.1^  A  schematic  diagram  of  the  programmable  matched 
filter  is  shown  in  Fig.  IV-4.  It  consists  of  a  300-stage  CCD  shift  register  with  parallel  outputs 
connected  to  an  array  of  300  sampling  fingers  on  a  p-type  silicon  substrate  which  is  in  close 
proximity  (typically  350  nm)  to  a  LiNbO^  delay  line.  The  RF  piezoelectric  fields  associated 
with  the  surface  wave  are  picked  up  by  the  sampling  fingers,  each  of  which  has  a  DC  bias  pre¬ 
viously  set  by  the  CCD-loaded  reference  pattern.  Nonlinearities  in  the  silicon  effectively  mul¬ 
tiply  the  RF  and  DC  potentials,  and  the  local  product  terms  are  capacitively  coupled  to  and 
summed  by  the  output  electrode,  thereby  producing  the  desired  cross-correlation  output. 
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Fig.  IV-4.  Schematic  diagram  of  a  SAW/CCD  programmable  matched  filter. 

This  output  electrode  is  a  0.25-mm-wide  A1  strip  deposited  on  the  LiNbOj  outside  the 
acoustic  beam  so  as  not  to  distort  the  SAW  phase  fronts.  Its  height  is  approximately  50  nm 
less  than  the  spacer-rail  height;  thus  it  is  gap-coupled  to  the  sampling  fingers.  Incorporation 
of  the  output  electrode  is  the  one  fabrication  detail  required  for  the  new  device  that  was  not 
needed  for  the  buffer  memory.  Electrical  differences  between  the  modes  of  operation  for  the 
devices,  however,  are  substantial. 

The  acoustoelectric  interaction  utilized  for  the  buffer  memory  is  a  linear  one,  while  that 
employed  in  the  matched  filter  is  nonlinear.  Furthermore,  as  discussed  below,  the  matched 
filter  is  a  slow-in,  fast-out  device,  in  contrast  to  the  fast-in,  slow-out  buffer  reported  earlier. 

The  sequence  for  matched -filter  operation  can  be  explained  by  reference  to  Fig.  IV-4. 

Prior  to  the  application  of  an  input  signal  the  sampling  fingers  are  precharged  to  an  appropriate 
potential  by  means  of  the  linear  array  of  enhancement-mode  MOS  transistors,  shown  at  the  top 
of  Fig.  IV-4.  After  the  biasing  transistors  are  turned  off,  allowing  the  sampling  fingers  to  float, 
a  reference  voltage  waveform  is  applied  across  the  input  gates  while  the  input  diode  is  pulsed 
at  the  CCD  clock  rate,  allowing  analog  samples  of  reference  charge  to  enter  the  CCD  in  the 
familiar  "fill-and-spill"  scheme.17  Serial  clocking  of  the  shift  register  is  halted  once  the 
300  <pi  wells  of  the  CCD  contain  the  desired  spatial  charge  pattern.  Charge  from  each  <f  well 
is  then  transferred  through  a  set  of  transfer  gates  to  the  well  associated  with  the  corresponding 
sampling  finger,  resulting  in  a  spatial  variation  of  the  initially  uniform  potential  on  the  array 
of  sampling  fingers. 

In  effect  this  process  biases  300  independent  varactors  which  perform  the  required  multi¬ 
plication  between  the  reference  charges  and  the  RF  fields  that  accompany  a  SAW  signal  as  it 
propagates  through  the  3.5-Msec-long  interaction  region.  The  silicon  nonlinearities,  previously 
exploited  in  the  SAW/CCD  accumulating  correlator,1 5  are  characteristic  of  both  the  drain  regions 
of  the  bias  transistors  as  well  as  the  MOS  regions  underlying  the  sampling  fingers.  Optimization 
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of  the  nonlinear  acoustoelectric  interaction  is  afforded  by  adjustment  of  the  uniform  finger  bias 
(as  described  earlier)  to  control  the  space-charge  layer  underlying  the  sampling  array.  The 
cross  correlation  between  the  SAW  signal  and  CCD  reference  is  obtained  by  summing  the  re¬ 
sultant  local  products  on  the  output  electrode.  For  the  special  case  in  which  the  SAW  signal  is 
an  impulse,  the  output  will  represent  a  time-compressed  serial  scan  of  the  reference  pattern, 
in  effect  producing  a  slow-in,  fast-out  buffer  memory. 

This  collinear  parametric  interaction  may  be  described  in  frequency-wavevector  (o>,  k) 

1 8 

space  simply  as  the  product  of  the  (a >,  k)  signal  wave  with  the  (0.  -k)  reference  pattern  to 
produce  a  correlation  at  (u>,  0).  The  matching  of  the  wavevectors  over  a  40-MHz  signal  band¬ 
width  is  achieved  with  precisely  the  same  sampling  scheme  as  was  implemented  and  reported 
for  the  buffer  memory.  In  brief,  the  CCD  cell  size  of  40.6  pm  defines  an  effective  sampling 
frequency  f  of  86.7  MHz  for  the  acoustoelectric  interaction,  and  therefore  the  signal  bandwidth 
which  can  be  match-filtered  without  sampling  ambiguities  is  in  the  range  85.7  to  128.6  MHz 
(fg  to  3/2  fg).  In  order  to  process  the  information  contained  on  a  SAW  input  at  an  RF  carrier 
f  within  this  band,  the  reference  pattern  on  the  300  fingers  must  be  at  the  matcing  spatial  fre¬ 
quency  (fc  -  f g ) / (3480  m/sec)  which,  for  the  100-kHz  clock  rate  used  in  the  CCD  shift  register, 
requires  an  audio  reference  waveform  in  the  0-  to  50-kHz  range. 

The  prototype  device  has  been  programmed  and  successfully  operated  as  a  matched  filter 
for  a  variety  of  signals  including  CW  tone-bursts,  linear-frequency-modulated  waveforms,  and 

biphase -modulated  Barker  codes.  Performance  of  the  prototype  device  to  data  has  been  most 
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adequately  tested  by  a  13-bit  Barker  sequence  (  +  +  +  +  +  --  +  +  -  +  -4)  which  would  produce,  in  an 

ideal  matched  filter,  a  symmetric  cross-correlation  envelope  consisting  of  a  strong  central 
peak  and  12  sidelobes  with  a  peak-sidelobe  amplitude  ratio  of  1/13.  The  actual  device  perfor¬ 
mance  is  shown  in  Fig.  IV-5,  which  is  an  oscillograph  of  the  cross  correlation  of  a  3.25-psec- 
long  13-bit  Barker  code  signal  on  a  107.2-MHz  carrier  frequency  with  a  13-bit  Barker  code 
reference  initially  loaded  into  the  CCD  on  the  matching  25-kHz  carrier.  For  the  signal  input 
power  at  the  1-dB  acoustic  compression  point  of  27  dBm,  the  correlation  peak  is  at  a  power 
level  of  -41  dBm.  which  is  substantially  above  the  —  98-dBm  thermal  noise  in  a  40-MHz  band¬ 
width.  In  future  devices  the  peak  output  power  could  be  increased  by  10  dB  and  the  dynamic 
range  raised  over  the  present  57-dB  limit  by  increasing  the  CCD  charge -handling  capability  so 


Fig.  IV-5.  Oscillograph  of  matched  filter 
operation  for  a  13-bit  Barker-encoded 
signal  3.25  psec  in  length.  Signal  carrier 
frequency  is  107  MHz  while  the  CCD- 
controlled  reference  waveform  was  en¬ 
tered  at  a  100-kHz  rate. 
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that  it  can  more  effectively  charge  the  sampling  fingers.  (Currently  the  finger-to-#>1  well 
capacitance  ratio  is  about  two.)  Although  the  3-dB  full  width  of  the  peak  matches  the  250-nsec 
bit  time  as  expected,  the  sidelobes  are  irregular  with  poorly  defined  nulls,  and  the  peak-to- 
sidelobe  ratio  is  only  -14.5  dB  as  compared  with  the  theoretical  -22.3  dB. 

Such  signal -dependent  spurious  sidelobes  are  not  uncommon  in  acoustoelectric  devices, 
and  are  typically  due  to  the  detection  by  the  silicon  circuit  of  bulk  acoustic  waves  gem  rated  at 
the  input  transducer.20  Spurious  sidelobes  due  to  such  bulk  waves  can  be  suppressed  to  a  level 
at  least  35  dB  below  the  correlation  peak  by  appropriate  design  of  the  delay  line.  Such  sup¬ 
pression  techniques  have  been  successfully  implemented  in  a  Schott ky -diode /LiNbO,  memory 
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correlator  structure  and  need  not  be  discussed  further  here. 

For  the  prototype  SAW/CCD  an  additional  spurious  signal  arises  from  a  transducer¬ 
like  detection  of  the  SAW  signal  by  the  sampling  finger  array.  This  response  is  predictably 
of  sinx/x  character,  where  x  =  300  jr(fc  -  f  )/f  is  the  normalized  frequency  excursion; 
for  fc  <  94  MHz  it  contributes  spurious  sidelobes  above  those  caused  by  bulk-wave  gener¬ 
ation.  Clearly,  an  optimized  matched-filter  device  would  require  suppression  of  this  non¬ 
programmable  transducer  effect  which,  near  the  sampling  frequency,  competes  with  the  CCD 
programmable  filter  output. 

Finally  there  is,  in  addition,  a  spurious  response  not  shown  in  the  oscillograph,  resulting 
from  direct  capacitive  coupling  between  the  signal  transducer  and  the  output  electrode.  In 
future  devices  this  unwanted  signal  can  be  reduced  to  an  insignificant  level  by  introducing  suf¬ 
ficient  physical  separation  between  the  input  transducer  and  the  output  electrode. 

Ultimately  the  quality  of  matched -filter  operation  will  be  dependent  upon  both  the  suppresion 
of  spurious  sidelobes  discussed  above  as  well  as  the  maintenance  of  phase  and  amplitude  fidelity 
of  the  reference  pattern  retained  by  the  sampling  fingers.  Holding  time  for  the  reference  is 
limited  by  the  thermal  generation  or  leakage  current  in  the  n+p  drain-substrate  junction  of  the 
bias  transistors.  Decay  of  the  reference  leads  to  a  3-dB  loss  in  the  correlation  peak  1  sec 
after  programming.  This  device,  therefore,  offers  an  excellent  duty  factor  with  only  6  msec 
currently  required  for  refreshing  or  reprogramming  the  reference.  Because  CCD  clocking 

rates  of  5  MHz  are  not  unreasonable,  the  writing  time  would  in  fact  be  reduced  to  about  120  usee. 

22 

As  discussed  elsewhere,  inhomogeneities  in  the  silicon  generate  fixed-pattern  noise  and  hence 
corrupt  the  reference  potentials.  A  quantitative  measure  of  the  small  degradation  in  signal 
processing  gain  due  to  such  defects  is  not  possible  in  the  prototype  device,  as  the  processing 
gain  is  limited  by  bulk-wave  effects. 

The  performance  of  this  SAW/CCD  programmable  matched  filter  is  given  in  Table  IV-1. 

In  summary,  we  have  demonstrated  an  important  analog  signal  processing  function  in  an  inte¬ 
grated  SAW/CCD  structure.  The  main  feature  of  this  device  is  that  it  can  correlate  wide  band¬ 
width  analog  signals  of  up  to  300  time -bandwidth  products  against  electronically  controlled 
reference  waveforms.  This  unique  hybrid  structure  makes  possible  real-time  processing  of 
wideband  inputs  combined  with  low-speed,  CCD  programmable  reference  patterns.  The  low- 
speed  CCD  reference  input  is  compatible  with  high-density,  low-cost  integrated  circuits  for 

storing  a  set  of  reference  waveforms.  _  „  , 

R.  W.  Ralston 

D.  L.  Smythe 
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TABLE  IV-1 

SAW/CCD  PROGRAMMABLE  MATCHED  FILTER 

Signal  Input  Center  Frequency 

107  MHz 

Signal  Input  Bandwidth 

40  MHz 

CCD  Clock  Rate 

100  kHz 

Reference  Input  Center  Frequency 

25  kHz 

Reference  Input  Bandwidth 

50  kHz 

Signal  Input  Duration 

3.5  p»ec 

Reference  Holding  Time 

1 . 0  sec 

Maximum  Input  Power 

27  dBm 

Maximum  Output  Power 

-41  dBm 

Dynamic  Range  over  kTB 

57  dB 

Spurious  (bulk  wave)  Sidelobes 

- 14  dB 

Time-Bandwidth  Product 

300 

C.  CCD  PROGRAMMABLE  TRANSVERSAL  FILTER 

A  16-tap  CCD  transversal  filter  (Fig.  IV-6)  with  tap  weights  programmable  as  4-bit  digital 
words  has  been  fabricated  and  tested.  The  device  has  been  successfully  used  as  a  programmable 
matched  filter  for  a  chirp  signal.  The  matched -filter  output  has  a  dynamic  range  of  50  dB. 

The  tap  weights  (h^)  in  this  device  are  represented  in  2's  complement  notation,  i.e., 

3 

hn=  2  ak>n2-k  —  a0,n  n  =  0,  1,  ...  15  (IV-1) 

k=l 

where  k  is  the  bit  position  (k  =  0  is  the  most  significant  bit  and  is  handled  separately  as  de¬ 
scribed  below),  n  is  the  tap  number,  and  a^  n  is  the  weighting  of  the  k**1  bit  of  the  nth  tap  word 
and  is  equal  to  0  orl.  The  filter  output  can  be  written  as 

15  3  15 

Vout<mTc>*  2  l  ak,n2'kvinf(m-n)Tc]-  l  aonVin  ><m  "  ">  Tc>  (,V-2' 

n=0  k=l  n=0 

where  T  =  l/f  ,  f  is  the  CCD  clock  frequency,  and  m  is  a  positive  integer. 

c  c  c  23  24 

The  filter  structure  is  based  on  a  device  described  previously  ’  and  is  shown  schemat¬ 
ically  in  Fig.  IV-7.  Each  reference  bit-word  is  serially  loaded  into  a  shift  register  and  stored 
in  the  latch  circuits.  The  output  of  each  triangular  section  represents  the  correlation  of  an 
analog  input  signal  with  a  specific  bit  of  the  reference  words.  As  can  be  seen  in  Fig.  IV-7,  the 
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Fig.  IV-6.  Photograph  of  4-bit,  16-tap 
CCD  programmable  transversal  filter. 
Each  triangular  section  is  one  bit. 


correct  filter  output  is  achieved  by  binary  weighted  summation  of  signals  from  three  sections. 
This  summation  can  be  performed  in  off-chip  analog  circuitry,  but  at  the  expense  of  additional 
components  and  power  dissipation.  A  novel  technique  for  performing  on-chip  weighted  signal 
summation  in  the  charge  domain  at  the  output  of  each  section  has  been  developed  and  is  de¬ 
scribed  here. 

The  operation  of  the  weighted  charge  summation  circuit  is  illustrated  in  Figs.  IV-8  and  -9. 
Each  triangular  section  has  a  long  charge  collection  node  along  the  outputs  of  the  delay  lines. 
This  collection  node  consists  of  two  diode  diffusions  (I>1  and  D2  in  Fig.  IV-8)  separated  by  a 
gate,  called  a  partition  gate,  forming  transistor  Qp.  The  operation  of  this  circuit  is  shown  in 
three  stages  beginning  [Fig.  IV-9(a) ]  at  a  time  when  charge  packets  in  the  last  <p2  storage  wells 
are  ready  to  be  transferred  to  the  collection  diodes  and  the  partition  gate  is  off.  When  the  <fi2 
clock  and  transfer  gate  go  low,  the  charge  from  the  16  delay  lines  is  transferred  to  the  two  col 
lection  diodes.  Since  each  collection  diode  services  a  different  set  of  delay  lines,  the  partition 
gate  must  go  on  to  allow  charge  exchange  and  to  equalize  the  potential  on  the  two  diodes 
[Fig.  IV-9(b) ].  The  partition  gate  then  turns  off,  the  transfer  gates  (Qtl  and  Q(2  in  Fig.  IV-8) 
go  on,  and  the  separated  quantities  of  charge  are  transferred,  one  to  the  charge  "dump"  and 
the  other  to  an  output  diode.  The  portion  Qg  [Fig.  IV-9(c) )  transferred  to  the  output  diode  is 
the  desired  signal,  and  is  related  to  the  total  charge  QT  from  the  triangle  section  by 
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Fig.  IV-7.  Schematic  describing  operation  of  the  programmable 
transversal  filter  shown  in  Fig.  IV-6. 
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Fig.  IV-8.  Schematic  of  the  on-chip,  binary-weighted  charge 
summation  scheme. 


CHARGE  COLLECTION  NODE 


PARTITION  GATE 

/  TRANSFER  GATE 


OUTPUT  DIODE 


Fig.  IV-9.  Cross  section  of  a  section 
of  the  CCD  programmable  filter  with  a 
corresponding  potential  diagram  dem¬ 
onstrating  the  charge  division  technique, 
(a)  Charge  packets  in  the  last  <p 2  wells 

are  ready  to  be  transferred  to  the  col¬ 
lection  diodes,  (b)  Partition  gate  is  on 
and  charge  is  equalized  on  the  two  col¬ 
lection  diodes,  (c)  Signal  charge  Qs  is 
transferred  to  the  output  diode. 
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where  C.  and  C,  are  the  capacitances  of  each  diode  segment.  The  capacitance  ratio  C,/(C .  +  C,) 

1  ^  4  12 

is  made  equal  to  2  in  accordance  with  Eq.  (IV-2).  The  charge  from  two  of  the  three  remaining 
sections  is  similarly  weighted  and  transferred  to  the  same  output  diode  where  it  is  sensed  by  a 
conventional  floating  diffusion  output  circuit  consisting  of  output  and  reset  transistors,  shown 
as  Qq  and  Qr  in  Fig.  IV-8.  In  the  present  structure  the  MSB  output  is  brought  out  separately 
to  be  subtracted  in  an  off-chip  differential  amplifier  from  the  output  of  the  remaining  three  sec¬ 
tions.  An  examination  of  Fig.  IV-6  shows  the  location  of  the  partition  along  the  staircase-shaped 
collection  diodes  at  the  outputs  of  all  but  the  MSB  triangle. 

To  demonstrate  programmable  matched  filtering  for  a  chirp  signal,  the  tap  weights  of  the 
device  have  been  progi  ammed  as  a  cosine  down-chirp,  i.e., 

hR  =  cos  8ir  t15^")2  n  =  0,  ...15 

The  impulse  response  of  the  device  is  shown  in  Fig.  IV-10(a).  When  a  negative  cosine  up-chirp 

•  V.  (m)  =  -cos  8ir  (77-)^  m  =  0,  ...  15 

in  10 

is  applied  to  the  CCD  signal  gate,  the  convolution  output  of  the  filter  has  the  correct  peak-to- 
sidelobe  ratio  for  this  waveform,  as  shown  in  Fig.  IV-10(b).  Operating  in  this  mode,  a  dynamic 

range  of  50  dB  has  been  measured.  .  . 

A.  M.  Chiang 

B.  E.  Burke 


Fig.  IV-10.  Programmable  matched  filtering 
using  the  4-bit,  16-tap  device,  (a)  Impulse 
response  of  the  device  with  tap  weights  pro¬ 
grammed  as  a  cosine  down  chirp,  (b)  Output 
of  the  device  when  a  cosine  up-chirp  is  ap¬ 
plied  to  inputs  of  the  CCD  and  the  tap  weights 
are  programmed  as  in  (a).  Output  is  the  con¬ 
volution  of  the  two  signals  and  gives  the  cor¬ 
rect  psak-to-sidelobe  ratio  for  this  waveform. 


Fig.  IV-11.  Diagram  of  system  for  measurement  of  CCD  quantum  efficiency. 
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Fig.  IV-12.  Details  of  gate  structure  of  the  CCD  showing  (a)  cross-sectional 
view  and  (b)  top  view  of  a  unit  imaging  cell. 


D.  CHARGE-COUPLED  DEVICES:  IMAGERS 


Measurements  of  the  spectral  quantum  efficiency  of  the  100-  x  400-element  CCD  imager 
have  been  made.  The  device  is  a  polysilicon  gate  front  illuminated  silicon  CCD  made  for  the 
GEODSS  (Ground  Electro -Optical  Deep  Space  Surveillance)  program  for  use  as  a  ground-based 
sensor  for  satellite  surveillance.  For  this  application  the  sensor  is  to  have  maximum  respon- 
sivity  to  an  air  mass  2  (AM2)  solar  spectrum  which  has  a  maximum  photon  flux  at  about  700  nm. 
Measured  data  span  the  wavelength  range  from  400  to  1100  nm,  and  typically  average  around 
50-percent  quantum  efficiency  over  a  midband  range  from  600  to  900  nm. 

The  quantum  efficiency  was  measured  by  illuminating  the  device  with  radiation  from  a 
monochromator  and  measuring  the  photocurrent  collected  by  the  CCD.  An  automated  system 
controlled  by  a  Hewlett-Packard  9825A  calculator  was  used  for  the  measurements  and  data 
reduction,  and  this  system  is  illustrated  by  the  diagram  in  Fig.  IV-11.  The  calculator  controls 
the  measurements  through  an  I/O  instrument  (Multiprogrammer)  containing  circuit  cards  to 
provide  the  interface  between  the  calculator  and  the  measurement  instruments.  A  card  provides 
the  timing  pulses  to  control  the  stepper  motor  which  sets  the  monochromator  wavelength.  The 
monochromator  employs  a  holographic  grating  which  is  usable  from  300  to  1100  nm.  The  slits 
used  in  these  measurements  produce  a  3-nm  optical  bandwidth.  A  portion  of  the  optical  beam 
is  reflected  by  the  beam  splitter  to  a  silicon  photovoltaic  detector.  The  detector  is  packaged 
with  an  op-amp  which  provides  a  voltage  output  proportional  to  the  optical  power  and  is  cali¬ 
brated  by  the  manufacturer.  The  photoelectrons  collected  by  the  CCD  are  transported  to  an 
output  sense  circuit  and  come  off  the  chip  at  a  diode  called  the  reset  diode.  The  photocurrent 
flowing  into  this  diode  is  measured  by  an  electrometer  from  which  a  voltage  proportional  to 
this  current  is  available.  The  detector  and  electrometer  voltages  are  multiplexed  through  re¬ 
lays  to  a  12-bit  A/D  converter  contained  in  the  Multiprogrammer.  The  optical  beam  image  is 
contained  entirely  within  the  active  areas  of  both  detector  and  CCD,  and  therefore  no  area  cor¬ 
rections  are  needed  in  the  data.  The  CCD  responsivity  is  thus  the  photocurrent  divided  by  the 
incident  optical  power.  The  CCDs  used  in  these  measurements  generally  had  room  temperature 
dark  currents  of  less  than  2  nA  (corresponding  to  6  nA/cm2),  while  saturation  level  photocur¬ 
rents  were  80  nA  or  more  at  the  clock  rates  used.  Because  of  this  high  signal -to -background 
ratio,  synchronous  detection  methods  were  not  required. 

The  CCD  structure  in  cross  section  is  shown  in  Fig.  IV-12(a),  and  a  top  view  of  a  unit  pixel 
in  Fig.  IV-12(b).  The  channel  stops  are  a  relatively  inactive  region  optically  but  account  for 
only  16  percent  of  the  area.  Thus  the  structure  has  a  high  percentage  of  active  area.  For  the 
device  discussed  below,  the  first-  and  second-level  polysilicon  thicknesses  were  3500  and  2600  A, 
respectively,  the  gate  SiC>2  and  Si3N4  thicknesses  were  both  800  A,  and  the  top  protective  oxide 
was  1.5  pm  thick.  The  first  polysilicon  was  deposited  in  an  RF-heated  reactor,  and  its  thickness 
is  somewhat  greater  than  desired  due  to  problems  of  uniformity  control  in  that  system.  A  low- 
pressure  hot-wall  CVD  (LPCVD)  deposition  system  became  available  for  the  second-level  poly¬ 
silicon.  This  system  has  superior  uniformity  and  thickness  control,  and  permitted  a  much 
thinner  second  polysilicon  layer  to  be  used.  This  is  an  important  consideration  for  the  CCD 
imager  because  the  polysilicon  thickness  has  a  major  influence  on  the  quantum  efficiency,  not 
only  in  the  blue  portion  of  the  spectrum  where  the  optical  attenuation  in  silicon  is  high,  but  also 
at  longer  wavelengths  where  optical  tailoring  techniques  can  be  used  to  minimize  reflection 
losses.25 
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A  representative  set  of  data  is  shown  in  Fig.  IV-13.  The  quantum  efficiency  averages 
approximately  50  percent  in  the  range  from  600  to  900  nm  for  this  device.  The  monochromator 
radiation  is  nearly  collimated  and  normal  to  CCD  surface.  In  an  optical  system  the  radiation 
will  be  incident  on  the  sensor  over  a  range  of  angles,  the  maximum  angle  from  normal  being 
arctan  f(l  +  d/D)/2f].  where  d  is  the  longest  dimension  of  the  sensor,  D  is  the  diameter,  and 
f  the  focal  length  of  the  optics.  In  the  GEODSS  systems  with  approximately  F/3  optics  this 
angle  will  be  about  10°.  Measured  data  show  very  little  variation  of  quantum  efficiency  up  to 
angles  of  15°.  The  variation  in  total  responsivity  and  average  quantum  efficiency  as  defined 
below  is  less  than  3  percent  over  this  range  of  angles. 

A  theoretical  calculation  of  the  quantum  efficiency  is  also  shown  in  Fig.  IV-13  based  on  the 

parameters  of  the  device.  The  calculation  uses  standard  methods  of  calculating  transmission 

26  27  28 

through  multiple  films,  and  the  optical  constants  of  single  crystal  silicon  '  are  used  for 

the  polysilicon.  The  formulation  of  the  carrier  generation  and  collection  is  based  on  the  work 
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of  Brown  and  Chamberlain.  The  most  striking  feature  in  the  comparison  of  experiment  and 

theory  is  the  relative  smoothness  of  the  data  compared  to  the  large  amplitude  peaks  and  valleys 

25 

of  the  theory.  Similar  results  have  been  noted  elsewhere  for  polysilicon  gate  imagers.  The 
discrepancy  has  been  attributed  to  the  effects  of  optical  scattering  at  the  polysilicon  grain 
boundaries. 

For  any  imaging  sensor  it  is  necessary  to  know  how  the  responsivity  of  the  sensor  matches 
the  radiation  output  of  the  target  or  image  being  viewed.  In  the  GEODSS  system  the  targets  of 
interest  are  satellites  which  are  viewed  in  reflected  sunlight.  Satellites  of  particular  interest 
are  those  in  geosynchronous  orbits,  and  it  is  assumed  that  the  spectral  irradiance  of  such  tar¬ 
gets  at  future  GEODSS  sites  will  be  approximately  that  of  an  AM2  solar  spectrum.  Efforts  to 
improve  the  CCD  quantum  efficiency  will  therefore  be  directed  toward  maximizing  the  respon¬ 
sivity  to  this  irradiance  distribution. 

A  plot  of  the  relative  AM2  solar  photon  flux  for  wavelengths  from  400  to  1100  nm  is  shown 
in  Fig.  IV-14.  The  peak  photon  flux  occurs  at  about  700  nm  and  has  an  overall  shape  similar  to 
the  CCD  quantum  efficiency.  Two  parameters  we  use  to  characterize  the  CCD  response  to  a 
particular  source  are  the  total  responsivity  R.rp  and  a  weighted  average  quantum  efficiency  rf 
They  are  defined  by  the  following  expressions: 
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where  R(A)  =  (qX/hc)  »j(A)  is  the  sensor  spectral  responsivity  in  A/W,  tj(X)  is  the  spectral  quantum 
efficiency,  and  h  =  Planck's  constant.  The  spectral  power  flux  W(A)  and  photon  flux  <p(A)  of  the 
source  are  related  by  <p( X)  =  (X/hc)  W(A).  The  average  quantum  efficiency  if  is  defined  for  the 
wavelength  range  of  400  to  1000  nm,  which  represents  the  useful  limits  of  the  CCD  sensitivity. 

Of  the  total  solar  AM2  energy,  approximately  70  percent  is  contained  within  this  band.  The 
CCD  response  beyond  1000  nm  is  not  considered  useful  due  to  radiation  from  night  sky  background. 
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For  the  data  of  Fig.  IV-13,  the  solar  AM2  responsivity  RT  is  153  mA/W  and  the  average 
quantum  efficiency  rj  is  41  percent,  while  the  theoretical  values  are  161  mA/W  and  41  percent, 
respectively.  This  close  agreement  is  somewhat  surprising  in  view  of  the  dissimilarities  of 
the  two  curves. 

Improvements  in  the  quantum  efficiency  can  be  expected  in  future  devices  due  to  the  im¬ 
proved  control  over  polysilicon  thickness  in  the  LPCVD  system.  Experiments  are  under  way 
to  make  devices  having  the  gate  structure  in  which  an  additional  silicon  nitride  layer  is  placed 
on  top  of  the  polysilicon2^  and  all  layer  thicknesses  adjusted  to  maximize  the  solar  AM2 
responsivity. 

J  B.  E.  Burke 
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V.  ANALOG  DEVICE  TECHNOLOGY 


A.  SUPPRESSION  OF  SPURIOUS  SIGNALS  IN  MEMORY  CORRELATORS 

The  development  of  an  acoustoelectric  memory-correlator  subsystem  has  been  reported.1’2 
The  heart  of  this  subsystem  is  a  memory  correlator  which  was  described  in  some  detail  in  an 
earlier  report.1  To  achieve  optimum  performance  of  the  subsystem,  and  any  system  of  which 
the  memory  correlator  is  a  part,  the  memory  correlator  must  be  as  free  of  spurious  signals  as 
possible.  Considerable  improvement  has  been  made  in  the  suppression  of  the  unwanted  spurious 
signals,  and  techniques  have  been  developed  which  promise  still  further  suppression. 

The  first  step  in  eliminating  spurious  signal  3  within  memory  correlators  was  to  identify  the 
source  of  these  signals  and  then  to  develop  suppression  techniques.  Four  separate  sources  have 
been  identified.  They  are:  electromagnetic  feedthrough  from  the  input  transducer,  bulk-wave 
generation  by  the  input  transducer,  bulk-wave  scattering  occurring  when  the  surface  wave  im¬ 
pinges  on  the  support  structure  for  the  silicon  strip,  and  edge  detection  by  the  edges  of  the 
silicon  strip.  The  problems  of  electromagnetic  feedthrough  and  bulk-wave  generation  by  the 
transducer  were  discussed  in  a  previous  report.1  New  approaches  have  been  taken  to  achieve 
suppression  of  the  other  spurious  signals. 

The  most  troublesome  of  the  spurious  signals  has  been  edge  detection  caused  by  the  stray 
capacitance  at  the  edges  of  the  silicon  strip.  The  result  is  an  output  signal  at  the  same  fre¬ 
quency  as  the  surface  wave  occurring  when  the  wave  first  passes  under  the  silicon  and  again 
when  it  leaves.  A  variety  of  methods,  including  slanting  the  ends  of  the  silicon  strip  and  using 
metal  shields  on  the  crystal  surface  under  the  ends  of  the  silicon  strip,  have  been  tried  with 
limited  success.  The  method  which  has  proven  successful  has  been  the  use  of  a  double - 
transducer  structure  in  conjunction  with  metal  shields,  as  shown  schematically  in  Fig.  V-i.  The 
double  transducer  launches  two  parallel  acoustic  beams  which,  because  of  the  180°  hybrid  at  the 
input,  are  out  of  phase  by  180*.  Identical  matching  networks  are  required  on  each  section  of  the 
double  transducer.  Ideally,  the  signals  detected  by  the  ends  of  the  silicon  would  be  out  of  phase 
and  cancel  completely.  In  practice  cancellation  is  incomplete,  but  the  use  of  metal  shields  is 
sufficient  to  suppress  any  residual  pickup.  The  use  of  the  double-transducer  configuration  and 
the  out-of-phase  drive  is  as  effective  in  suppressing  electromagnetic  feedthrough  as  the  bal¬ 
anced  device  which  was  described  earlier.1 

With  the  double  transducer  it  is  unwieldy  to  use  a  multistrip  coupler  to  eliminate  bulk  waves 
generated  by  the  transducers.  In  practice  it  has  been  observed  that  the  level  of  spurious  bulk- 
wave  signals  generated  by  the  double-transducer  configuration  is  low,  perhaps  because  any  bulk 
waves  generated  by  the  antiphase  drive  in  the  two  transducers  tend  to  cancel.  The  addition  of  a 
transparent  urethane  absorber  on  the  back  of  the  surface-wave  substrate  has  proven  capable  of 
eliminating  most  of  the  remaining  spurious  bulk  waves  generated  by  the  transducers.  Experi¬ 
ments  are  currently  in  progress  to  determine  if  the  127.86°  rotated  Y-cut  of  LiNbOj  will  improve 
this  situation.  This  cut  is  known  to  have  very  low  coupling  to  bulk  waves. 

The  last  remaining  source  of  spurious  signals  is  the  scattering  into  bulk  waves  which  occurs 
when  the  surface  waves  are  incident  on  the  rails  that  support  the  silicon  strip;  however,  this 
spurious  signal  has  been  greatly  reduced  by  the  use  of  the  double -transducer  configuration.  The 
surface  wave  launched  by  this  configuration  has  a  null  at  the  center  of  the  beam;  thus  scattering 
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Fig.  V-i.  Schematic  diagram  of  a  memory  correlator  containing  a  double  input  transducer 
and  metal  patterns  for  suppression  of  edge  detection. 
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Fig.  V-2.  Major  parts  of  a  memory  cor¬ 
relator  containing  double  input  trans¬ 
ducers.  Top:  Silicon  strip  mounted  on 
Kapton.  Middle:  LiNbC>3  substrates  and 
matching  circuits  mounted  in  package 
base.  Bottom:  LED  array  for  erasure. 


off  the  central  supporting  rail  is  eliminated.  However,  some  surface-wave  energy  is  incident 
on  the  outer  rails,  most  likely  as  the  result  of  diffraction  from  the  transducers,  and  hence 
some  small  amount  of  energy  is  still  converted  to  bulk  waves. 

A  preassembly  view  of  a  memory  correlator  with  the  double -transducer  configuration  is 
shown  in  Fig.  V-2.  The  LiNbOj  crystal  is  mounted  in  the  package  with  matching  network  at 
either  end.  The  double  transducer  is  visible  at  the  left-hand  end  of  the  crystal,  while  the  metal 
shields  that  are  located  under  the  ends  of  the  silicon  are  the  very  small  patterns  in  the  central 
region  of  the  crystal.  The  dark  band  that  covers  the  central  section  of  the  crystal  is  the  trans¬ 
parent  ground  plane  which  is  on  the  backside  of  the  LiNbO,.  Above  the  crystal  is  located  the 
silicon  strip  containing  the  Schottky  diode  array.  The  object  located  below  the  crystal  in 
Fig.  V-2  is  the  bottom  cover  plate,  which  contains  the  LED  array  that  is  used  for  erasure  of  the 
stored  pattern  in  the  diode  array. 

Figure  V-3(a)  shows  the  response  of  the  memory  correlator  when  used  to  process  a  5-psec- 
long  50-MHz-bandwidth  linear-FM  waveform.  The  approximately  (sinx)/x  sidelobes  due  to  the 
compression  of  an  unweighted  chirp  are  readily  app&i  »nt.  In  Fig.  V-3(b)  the  remaining  residual 
spurious  signals,  virtually  all  caused  by  bulk-wave  scattering,  can  be  observed.  (Note  the  in¬ 
creased  gain  of  20  dB. )  At  the  present  time  the  spurious  signals  introduce  sidelobes  which  are 


Fig.V-3.  Response  of  memory  correlator,  (a)  Compressed  pulse  obtained 
during  matched  filtering  of  a  5-psec  50-MHz  linear-FM  waveform,  (b)  Gain 
increased  by  20-dB  to  show  residual  spurious  signals  which  occur  at  the 
output  even  when  no  reference  is  stored. 


on  the  order  of  35-dB  below  the  peak  memory  correlator  output.  Experiments  are  currently 
under  way  in  an  attempt  to  eliminate  the  remnants  of  this  spurious  signal  by  using  a  low-bulk- 
wave  cut  of  LiNbQj  and  by  forming  a  guiding  structure  on  the  surface  of  the  LiNbOj  so  that  the 
surface-wave  energy  will  not  be  incident  on  the  supporting  rails. 

D.  H.  Hurlburt 


B.  MNOS  NONVOLATILE  ANALOG  MEMORY 

MNOS  capacitor  structures  are  being  investigated  for  use  as  nonvolatile  analog  memories. 
Such  a  memory  is  compatible  with  present  CCD  and  SAW  /CCD  technology  and  would  greatly 
extend  the  signal  storage  time  of  these  devices. 

In  the  past  few  years  a  widespread  technological  effort  has  been  directed  at  the  development 
of  MNOS  structures  for  the  storage  of  digital  data.  The  device  is,  in  its  simplest  form,  a  dual¬ 
dielectric  capacitor  cell  as  shown  in  Fig.  V-4.  It  consists  of  a  metal  gate  electrode  (Cr  in  our 


Fig.  V-4.  MNOS  charging  with  majority 
carriers.  Oxide  field  Eo  is  linearly 
dependent  on  the  applied  voltage  and  the 
quantity  of  stored  charge. 


case,  but  more  commonly,  Al),  a  chemically-vapor-deposited  silicon-nitride  layer  about  500  A 
thick,  a  thermally  grown  silicon-dioxide  layer  only  25  to  40  A  thick,  and  a  silicon  substrate 
(shown  as  n-type).  Consider  a  positive  voltage  applied  to  the  gate.  The  positive  bias  serves  to 
accumulate  majority  carriers  (electrons)  at  the  Si-SiOz  interface,  and  thus  there  is  little  poten¬ 
tial  drop  across  the  silicon.  For  moderate  gate  voltages  of  the  order  of  30  V,  there  exists  a 
sufficiently  intense  field  in  the  oxide  to  cause  tunneling  of  electrons  through  the  oxide.  The 
amorphous  nitride  layer  is  grown  with  very  low  conductivity,  and  hence  this  layer  serves  to 
block  most  of  the  electron  flow.  As  a  result,  the  electrons  are  captured  at  or  near  the  oxide - 
nitride  interface  by  a  high  density  of  traps  which  exist  at  the  interface  or  in  the  nitride.  The 
precise  nature  of  the  tunneling  current  and  trapped  charge  is  determined  by  the  spatial  and 
energy  distribution  of  empty  and  occupied  traps. 

Several  charge  transport  mechanisms  play  a  role  in  the  charging  and  discharging  of  these 
devices.  A  brief  examination  of  these  mechanisms  is  fruitful. 

The  band  diagram  for  majority  carrier  injection  is  shown  in  Fig.  V-5.  Two  tunneling  mech¬ 
anisms  are  most  plausible:  direct  tunneling  of  valence-band  electrons  to  traps  (Fig.  V-5(a)) 
and  band-to-band  (modified  Fowler-Nordheim)  tunneling  from  the  silicon  conduction  band  to  the 
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Fig.V-5.  Band  diagram  of  MNOS 
interfacial  region  during  majority 
charging. 


nitride  conduction  band  followed  by  trapping  |  Fig.  V-5(b)J .  In  both  cases,  the  tunneling  current 
is  exponentially  dependent  on  the  field  strength  and  oxide  thickness.  For  instance,  the  band- 
to-band  tunneling  current  may  be  written3 


JBB  =  CFNEoxPoxPN 


(V-l) 


where  Cp^  is  a  constant,  characteristic  of  Fowler-Nordheim  tunneling  from  silicon  into  Si02, 
and  Eqx  is  the  oxide  field  strength.  Pqx  and  PN  are  the  tunneling  probabilities  through  the 
oxide  and  nitride,  respectively,  given  by 

P  =  exp  {—  tT-  /  2qm*  [<p3/2  -  (<p ,  -  d  E  )3/2  ]/E  }  (V-2) 

ox  3h  si  ^  ox  lvl  ox  ox  J/  oxJ 

and 

PN  =  exp{~£;JZ*mN  {(pi  -^2-doxEox)V2/EN^  (V_3) 


where  <p ^  =  3.1  eV  and  <p2  =  1.05  eV  are  the  barrier  heights  (defined  in  Fig.  V-5),  dQx  is  the 

oxide  thickness,  and  m*  and  m*  are  the  electron  effective  masses  in  the  oxide  and  nitride, 
ox  N 

The  nitride  field  E^  is,  in  this  approximation,  assumed  to  be  related  to  the  oxide  field  by 
Ejj  =  (e  A™)  E  ;  i.e.,  the  space  charge  of  the  trapped  carriers  is  ignored. 

il  c/A  IN  DA 

Direct  tunneling  occurs  to  nitride  and/or  interfacial  traps  which,  in  the  absence  of  an  applied 
voltage,  are  located  energetically  in  the  nitride  bandgap  directly  opposite  the  silicon  gap;  other¬ 
wise  the  trapped  charges  would  quickly  back-tunnel  to  the  silicon.  Direct  tunneling  is  thought 

4  ° 

to  be  more  important  in  thin-oxide  devices  (425  A),  and  less  so  in  the  thicker-oxide  devices 
under  consideration. 

Carriers,  once  trapped,  may  escape  by  two  routes:  back-tunneling  through  the  oxide  to  the 
silicon,  or  conduction  (by  several  mechanisms)  through  the  nitride  to  the  gate  electrode.  The 
back-tunneling  current5  is  exponentially  dependent  on  the  oxide  thickness,  trap  location  in  the 
nitride,  trap  depth,  and  quantity  of  trapped  charge  (through  its  space-charge  effect). 

The  nitride  conduction  current,  on  the  other  hand,  is  critically  dependent  on  the  nitride 
processing  and  ambient  temperature.  Three  components  of  the  nitride  current  can  be  identified^1 
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Ji  =  cien  exP^  K  -(qEN  Acd)1/21AT}  (V.4) 

J2=C2ENexPf-E2/ENl  (V.5) 

J3  =C3ENexpf-q^3/kT]  (V.6) 

The  first  component  is  due  to  field -enhanced  thermal  excitation  of  trapped  electrons  into  the 
conduction  band  (Poole-Frenkel  effect).  The  constant  C.  and  dynamic  dielectric  constant 
are  dependent  on  the  nitride  deposition  and  processing  parameters.7' 8  The  second  component 
is  a  tunneling  current  resulting  from  field  ionization  of  trapped  electrons.  The  third  compo¬ 
nent  is  an  ohmic  current,  ascribed  to  the  hopping  of  thermally  excited  carriers  from  one 
isolated  state  to  another. 

8  At  r°°m  temPerature  the  component  dominates  for  electric  field  strengths  E  >  3  x 
10  V/m,  and  the  ohmic  component  dominates  for  weaker  fields,6  although  recent  experimental 
evidence  indicates  that  J4  dominates  even  for  E^  >  108  V/m  (Ref.  9).  The  component  J,  is 
important  only  at  low  temperatures.  Thus,  for  devices  operating  at  room  temperature  the 
Poole-Frenkel  component  may  be  expected  to  limit  the  value  of  the  charge  stored  during  a 
majority-carrier  writing  period  (when  nitride  fields  are  typically  6  X  108  V/m),  a  steady-state 
value  of  stored  charge  being  reached  when  the  nitride  current  equals  the  oxide  current.  This 
component  may  also  contribute  to  the  initial  loss  of  stored  charge;  however,  nitride  fields  in 
this  case  are  due  to  the  stored  charge  alone  and  are  typically  3  x  108  V/m  or  less.  The  ohmic 
current  may  be  expected  to  play  a  role  in  this  discharge. 

The  key  to  an  analog  memory  is  to  transfer  a  packet  of  minority  carriers  from  the  silicon 
to  the  nitride.  In  digital  devices,  minority  carrier  injection  is  achieved  by  providing  an  essen- 
lally  unlimited  supply  of  minority  carriers  to  an  inversion  layer,  by  means  of  a  field-effect- 
transistor  structure  or  avalanche  breakdown  in  the  depletion  layer.  Tunneling  through  the  oxide 
to  the  nitride  traps  then  proceeds  in  much  the  same  manner  as  for  majority-carrier  injection 
although  the  barrier  heights  and  trap  levels  (and  thus  time  constants)  may  be  different  For  ’ 
analog  memory,  a  measured  amount  of  minority  carriers  forms  the  inversion  layer.  Ultimately 
these  packets  of  carriers  will  be  clocked  under  the  memory  gate  in  CCD  fashion;  in  our  experi-  ' 
mental  devices,  the  earners  are  provided  by  injecting  photons  into  the  depletion  layer 

Minority-carrier  injection  is  illustrated  in  Fig.  V-6.  A  negative  gate  voltage  forms  a  de¬ 
pletion  layer  in  the  n-type  silicon.  In  the  absence  of  an  inversion  layer,  almost  all  the  applied 


Fig.  V-6.  MNOS  charging  with  minority 
carriers.  Oxide  field  E0  is  linearly 
dependent  on  inversion  layer  Charge  and 
only  weakly  dependent  on  applied  volt¬ 
age  and  stored  charge. 
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potential  (and  potential  of  any  trapped  charge)  is  dropped  across  the  depletion  layer;  the  oxide 
field  is  very  low.  The  injection  of  minority  carriers  collapses  the  depletion  layer  somewhat 
and  produces  an  oxide  field  which  is  nearly  proportional  to  the  quantity  of  minority  carriers. 
Ideally,  all  the  minority  carriers  would  tunnel  through  the  oxide  but,  because  the  tunneling  de¬ 
pends  on  the  self -field  of  the  minority  carriers,  a  certain  quantity  of  minority  charge  must 
remain  in  the  silicon.  This  necessitates  the  use  of  a  "fat  zero"  but  does  not  affect  the  linearity 
of  the  analog  memory. 

The  potential  barrier  is  of  the  same  form  as  given  by  Eqs.  (V-2)  and  (V-3),  but  the  pre¬ 
exponential  factor  in  Eq.  (V-l)  is  modified  to  account  for  the  finite  number  of  available  carriers. 
Additional  charge  transfer  is  effected  by  the  back-tunneling  of  stored  electrons,  induced  by  the 
increased  oxide  field. 

The  analog  storage  cell  is  operated  as  follows:  First,  a  positive  pulse  (=»30  V  and  some 
milliseconds  in  duration)  is  applied  to  the  gate.  This  causes  electron  tunneling  to  the  nitride 
traps  and  brings  the  device  to  a  known  "reset"  state.  A  negative  gate  pulse  (=<  —20  V)  is  then 
applied  to  deplete  the  silicon,  the  minority-carrier  packet  is  injected,  and  the  tunneling  of  mi¬ 
nority  charge  takes  place.  The  stored  charge  may  be  read  nondestructively  by  a  simple  C-V 
measurement;  the  flat-band  voltage  shift  is  proportional  to  the  stored  charge.  Ultimately,  in 
an  integrated  CCD-MNOS  device,  the  stored  signal  charge  would  be  read  nondestructively  by 
filling  the  depleted  potential  well  in  the  silicon  below  the  MNOS  gate. 

Here  we  wish  to  describe  the  encouraging  results  obtained  on  isolated  MNOS  memory  cells 
which  were  fabricated  on  10-fi-cm  n-type  Si  without  peripheral  input/output  circuitry.  The 
geometry  of  these  stand-alone  capacitors  is  depicted  in  Figs.  V-4  and  -6.  As  was  previously 
described,  a  positive  gate  voltage  is  first  applied  to  reset  the  memory  traps  to  a  large  electron 
occupation.  Then  a  comparable  negative  gate  voltage  is  applied  which  deep  depletes  the  silicon. 
Held  in  the  dark,  the  thermal  generation  of  minority  carriers  occurs  very  slowly  (full  inversion 
takes  2  min.),  and  thus  for  short  times,  insufficient  holes  and  insufficient  oxide  field  are  present 
for  appreciable  tunneling  to  occur.  Into  this  empty  deep-depleted  well  we  are  able  to  inject  a 
controlled  light  flux  by  illuminating  the  semi-transparent  Cr  gate  with  an  array  of  LEDs.  The 
photon  flux  absorbed  in  the  depletion  layer  generates  electron-hole  pairs;  the  depletion  field 
sweeps  the  electrons  into  the  neutral  bulk,  while  the  holes  are  swept  to  the  Si-SiO.,  interface 
and  there  form  an  inversion  layer.  The  formation  of  the  inversion  layer  causes  a  collapse  in 
the  depletion  field,  hence  a  redistribution  of  the  applied  voltage  occurs,  and  a  significant  tunnel¬ 
ing  field  is  obtained  in  the  oxide.  Holes  from  the  controlled  inversion  layer  are  able  to  occupy 
the  nitride  traps.  The  results  below  show  that  the  trapped  charge  is  linearly  proportional  to  the 
injected  inversion  charge  over  a  wide  range.  This  provides  a  suitable  analog  store  mechanism. 

The  results  of  our  experiment  are  shown  in  Fig.  V-7(a),  which  is  a  family  of  capacitance- 
voltage  characteristics  for  a  typical  MNOS  cell.  For  each  curve,  the  read  waveform  is  ramped 
from  +15  to  -15  V  in  100  ^sec  [right -to -left  in  Fig.  V-7(a)J,  which  sweeps  the  silicon  surface 
from  accumulation  through  flat-band  and  into  deep  depletion  with  a  corresponding  swing  in  ca¬ 
pacitance  from  about  1250  pF  to  nearly  zero.  The  read  waveform,  of  course,  is  not  sufficient 
to  produce  tunneling,  hence  it  does  not  perturb  the  trapped  memory  charge.  The  right-most 
curve  is  the  result  of  a  reset  operation  at  +  32  V.  All  other  curves  are  for  successively  increas¬ 
ing  amounts  of  signal  (inversion)  charge,  each  injected  optically  within  500  msec  at  -27  V  with 
intervening  resets. 


Fig.V-7.  (a)  Capacitance-voltage  characteristic  of  typical  MNOS  cell  indicating 
nondestructive  readout  of  different  analog  signal  charges,  (b)  Reduced  data  ob¬ 
tained  from  (a)  indicating  a  large  linear  range. 


Figure  V-7(b)  is  a  plot  in  which  we  indicate  the  midpoint  voltage  for  each  of  the  above  C(V) 
curves  as  a  function  of  integrated  light  flux  or,  equivalently,  the  analog  signal  charge.  Note 
that  over  a  voltage  swing  of  nearly  13  V  the  stored  charge  is  essentially  a  linear  (to  within  1  dB) 
replica  of  the  input  signal.  Certainly,  then,  large  quantities  of  charge  can  be  stored  in  a  linear 
analog  fashion.  This  maximum  linear  analog  charge-handling  capacity  we  label  as  AQm. 


VOLTAGE 

(c) 


Fig.V-7.  (c)  Capacitance-voltage  characteristic  around  the  central  point  in  (b)  indicating 
the  sensitivity  of  readout  to  small  increments  in  stored  charge. 


Figure  V-7(c)  is  the  expanded  region  about  the  central  point.  Here  the  smallest  voltage 
shifts  ( — 80  mV)  result  from  increments  in  signal  charge  equal  to  AQj^/220.  These  shifts  are 
clearly  discernible,  are  well  above  the  noise  level  in  this  read-out  scheme,  and  indicate  an 
analog  dynamic  range  for  this  single  MNOS  cell  in  excess  of  47  dB.  The  results  confirmed  our 
expectation  that  it  should  be  possible  to  enter  analog  signals  in  the  form  of  minority  carriers 
into  MNOS  storage. 

Of  major  concern  is  the  loss  of  signal  charge  with  time,  a  process  which  limits  the  useful 
storage  time  of  these  devices.  Figure  V-8  is  a  plot  of  experimental  measurements  of  this 
decay  in  our  devices.  The  decay  of  stored  charge  is  clearly  logarithmic  in  time,  although  for 
positive  flat-band  voltages  (negative  charge),  the  onset  of  this  decay  is  somewhat  delayed. 

(This  is  believed  to  be  due  to  the  low  oxide  field  existing  before  the  buildup  of  a  thermally  gen¬ 
erated  inversion  layer.)  The  rate  of  decay  (volts  per  decade)  is  evidently  dependent  on  the  ini¬ 
tial  flat-band  voltage.  A  plot  of  the  rate  of  decay  vs  initial  flat-band  voltage  (Fig.  V-9),  in  fact, 
reveals  a  linear  dependence,  as  has  been  predicted  theoretically.^  Thus,  although  the  decay  is 
logarithmic  in  time,  at  any  one  moment  a  linear  relation  holds  between  the  initial  analog  signal 
and  the  remaining  stored  charge.  The  magnitude  of  the  decay  rate  is  not  as  critical  as  its 
linear  dependence  on  initial  charge. 

Numerical  models  of  these  MNOS  capacitors  are  being  made  with  the  intention  of  gaining  a 
better  understanding  of  the  important  mechanisms  in  the  charging  and  discharging  of  these 
devices.  Further  improvements  in  the  charge  retentivity  of  the  memory  cells  are  also  depen- 

O  O 

dent  on  the  ability  to  grow  reproducible  thin  oxides  ( —  30  A,  with  an  error  of  2  A  or  less)  and  to 
deposit  low-conductivity  nitrides  with  reproducible  trap  levels  and  densities. 

In  summary,  we  have  demonstrated  nonvolatile  storage  of  analog  signal  levels  in  isolated 
solid  state  devices  with  a  dynamic  range  of  greater  than  40  dB.  Initial  results  indicate  that  the 
decay  rate  (volts  per  decade  of  storage  time)  is  linearly  dependent  on  the  initial  signal  level; 
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V„.  tvol»> 


Fig.  V-8.  Decay  of  stored  charge  (as  represented  by  the  midpoint 
voltage  VMp)  with  time  after  writing. 


Fig.V-9.  Hate  of  decay  of  stored 
charge  (midpoint  voltage),  in  units 
of  volts  per  decade  of  storage  tipne, 
as  a  function  of  initial  midpoint 
voltage.  Data  from  Fig.  V-8. 
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thus,  a  waveform  sampled  at  a  series  of  storage  sites  may  be  attenuated  in  time  but  not  dis¬ 
torted.  Even  so,  a  lower  decay  rate  would  be  advantageous.  Efforts  are  under  way  to  obtain 
a  better  understanding  of  storage  and  conduction  mechanisms  in  the  dual  dielectric  films.  Also, 
better  control  of  the  oxide  growth  and  nitride  deposition  processes  is  essential  to  improved 

charge  retention.  r.  w.  Ralston 

R.  S.  Withers 
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